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ABSTRACT 

The Fundamental Plane of early type galaxies is observed to have finite thickness and to be tilted 
from the virial relation. Both of these represent departures from the simple assumption that dynamical 
mass-to-light ratios (Mdyn/L) are constant for all early type galaxies. We use a sample of 16,000 
quiescent galaxies from the Sloan Digital Sky Survey to map out the variations in Mdyn/L throughout 
the 3D Fundamental Plane space defined by velocity dispersion (a), effective radius (Re), and effective 
surface brightness (/e). Dividing Mdyn/L into multiple components allows us to separately consider 
the contribution to the observed Mdyn/L variation due to stellar population effects, IMF variations, 
and variations in the dark matter fraction within one Re- Along the FP, we find that the stellar 
population contribution given some constant IMF (M^_/j\/i?/L) scales with a such that M-^jmf/L oc 
f{cT). Meanwhile, the dark matter and/or IMF contribution {Mdyn/Mi^jMp) scales with Mdyn such 
that Mdyn/Mi,jMF oc g{Mdyn)- This means that the two contributions to the tilt of the FP rotate 
the plane around different axes in the 3D space. The observed tilt of the FP requires contributions 
from both, with dark matter/IMF variations likely comprising the dominant contribution. Looking 
at Mdyn/L variations through the thickness of the FP, we find that Mdyn/L variations must be 
dominated either by IMF variations or by real differences in dark matter fraction with Rg. This 
means that the finite thickness of the FP is due to variations in the stellar mass surface density 
within Re (E^^ja/f), not the fading of passive stellar populations and it therefore represents genuine 
structural differences between early type galaxies. These structural variations arc correlated with 
galaxy star formation histories such that galaxies with higher Mdyn/Mi^jMF have higher [Mg/Fe], 
lower metallicities, and older mean stellar ages. We discuss several physical mechanisms that might 
explain the observed co- variation between Mdyn/Mi,jMF and galaxy star formation histories. It is 
difficult to explain the observed enhancement of a-elements in lower-surface-brightness galaxies by 
allowing the IMF to vary. Differences in dark matter fraction can be produced by variations in the 
"conversion efficiency" of baryons into stars or by the redistribution of stars and dark matter through 
dissipational merging. The former explanation, specifically a model in which some galaxies experience 
low conversion efficiencies due to premature truncation of star formation, provides a more natural 
explanation for the co- variation of Mdyn/M^^jMF and the observed stellar population properties. 

Subject headings: galaxies: elliptical and lenticular, galaxies: structure, galaxies: evolution 



1. INTRODUCTION 

Early type galaxies are observed to obey many scaling 
relations among their structural properties. Early work 
identified a number of ID relations, such as the Faber- 
Jackson relation bet ween galaxy luminosity (L) and ve- 
locity dispersion (a) (iFaber fc Jacksonlll976f ). variations 
in ga laxy mass-to-light ratio (M/L) versus L (jTinslevI 
119811 : [Fabcr ct_al. 1983j, and various correlations of a 
and L with galaxy eff ective radius (R, e.) and effective 
surface brightness (fie) ()Kormendvlll985f ). or with galaxy 
core r adius (re) and central surface brightness (/q) (jLaueii 
[19851) . 

These relations generally refiect a ID mass sequence 
of galaxies. It is clear however that early type galax- 
ies comprise at least a two-parameter family in terms 
of their structure. In the three-dimensional parameter 
space of galaxy properties defined by a, Rg, and surface 
brightness (expressed as fj,e in magnitudes or as Ig in lin- 



ear flux units), early type galaxies populate a relatively 
tight two-dime nsional plane, known as the Fundamen- 



tal P lane (FP, iDiorgovski fc Davis! [I987t [Dressier et all 
|1987| ). Projections of the FP appear narrow in some ori- 
entations, leading to the seemingly ID relations listed 
above. 

The FP can be understood as a manifestation of the 
virial plane predicted for relaxed systems, under the as- 
sumption that galaxy mass-to- light ratios (Mdyn/L) are 
constant (or at least smoothly varying) for all galaxies. 
If Mdyn/L were strictly constant and if structures were 
homologous for all early type galaxies, the FP would be 
equivalent to the virial plane, which takes the form Rf, ex 
a"^!'^, and would be infinitely thi n. Instead, the FP is 
rotated from the vi rial pl a ne (e.g., [Dressier et al.l 119871 : 
IDiorgovski fc David [19871) . lJ0rgensen et al.l ([1996[ ) find 



Re on a 



1.24 r-0.82 



for the FP of local cluster early type 
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galaxies, as measured in the Gun n r band. In addition , 
the FP has finite thickness (e.g., 'J orgensen et al.l 119961: 
Forbes et al.l[l99l fW uvts ct al. 20041: JHvde fc Bernardil 



20091 : iGargiulo et al.ll2009l) indicating further departures 
from a single virial plane. The "tilt" of the FP and its 
finite thickness indicate departures from the simplistic 
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assumption that Mdyn/L is constant for all early type 
galaxies and that their structures are homologous. 

This is the third in a series of four papers that ex- 
plores the mapping between the 2D family of early type 
galaxy stellar populations and their structural proper- 
ties. The first three papers in the series systematically 
gather together a number of galaxy properties that vary 
throug h the differen t dimen sions of galaxy parameter 
space. [Graves et al.l (|2009al . hereafter Paper I) showed 
that stellar populations vary differently as functions of 
a than as functions of L. This systematic variation sug- 
gests an und e rlying multi-dimensional parameter space. 
[Graves et al.l ()2009bl hereafter Paper II) mapped out the 
second dimension of stellar population properties in 3D 
Fundamental Plane space and correlated stellar popula- 
tion variations with FP parameters. In this paper, we 
focus on the corresponding mass-to-light ratios and how 
they relate to the distribution of galaxies in FP space. 

We show here that the stellar population variations 
measured in Paper II do not provide enough variation in 
galaxy mass-to-light ratios to explain either the tilt or 
the thickness of the FP. Some correlated variation is also 
required in either the initial mass function (IMF) with 
which galaxies form their stars, or in the dark matter 
fraction inside Re- At the same time, these variations are 
a key component of structural differentiation in galaxies 
that must also be related to their star formation histo- 
ries. The results of these first three papers are gathered 
together in Paper IV (Graves et al. 2010, submitted to 
ApJ), which uses them to propose scenarios for galaxy 
formation and evolution through the thickness of the FP. 

Section [2] briefly outlines the sample of galaxies used in 
this work, as well as our methods for grouping together 
similar galaxies to produce high signal-to-noise spectra 
that span the space of galaxy properties. In section[3l we 
discuss the various physical mechanisms that may con- 
tribute to the observed variations in Mdyn/ L through- 
out the FP. Some of this variation is due to stellar pop- 
ulation effects, but quantifying these depends critically 
on the stellar population models used to measure stellar 
mass-to-light ratios [Mi^jMF /L). We therefore investi- 
gate multiple ways of measuring M-^jmf/L from stellar 
population models in sectional This comparison allows 
us to explore the range of possible M^^jme/L values for 
our sample galaxies and to understand the biases inher- 
ent in each model. We demonstrate that, regardless of 
which method is chosen to measure M^^jme/L, the ob- 
served stellar population variations are inadequate to ex- 
plain either the tilt or the thickness of the FP. Section \5\ 
shows that the two contributions to the FP tilt (1: from 
known stellar population variations and 2: from IMF 
and/or dark matter fraction variations within i?e) rotate 
the plane around different axes in 3D Fundamental Plane 
space. It also investigates the various contributions to 
the thickness of the FP and shows that the spread in 
Mdyn/ L through the plane is due primarily to the sec- 
ond factor: variations in the IMF and/or dark matter 
fractions. Section [6] discusses physical models that may 
produce the observed correlations between dark matter 
fraction or IMF variation and galaxy star formation his- 
tories. It illustrates that the observed trends can be ex- 
plained by a model in which some galaxies experience low 
"conversion efficiency" turning baryons into stars due to 
trucated star formation histories. Finally, section [7] sum- 



marizes our conclusions. 

2. DATA 

The data used in this analysis are the same as those 
used in Papers I and II. They consist of a sample 
of ^^16,000 non-star-forming (quiescent) galaxies cho- 
sen from the SDSS spectroscopic Main Galaxy Sample 
(jStrauss et al.lf2002l ) in a relatively narrow redshift range. 
The sample selection is described in detail in Paper I. 

Briefly, we construct a sample of quiescent galax- 
ies by requiring that their spectra contain no de- 
tectable emission in either Ha or the [Oii]A3727 dou- 
blet (all measured line strengths are below a 2a detec- 
tion thres holdo, ba s ed on the emission line measure- 
ments o£^^^iU200i)_Jor SDSS Data Release 4 
(DR4. IAdelman-McGarthv et al.|[2006l) galaxies. We fur- 
ther limit the sample to a relatively narrow redshift slice 
(0.04 < z < 0.08) and to galaxies with light profiles 
consistent with early type morphologies. These criteria 
produce a sample of galaxies that fall on the red sequence 
in a color- magnitude diagram (see Figure 1 of Paper I), 
and which typically have bulge-dominated morphologies. 

The various FP parameters are obtained for all galax- 
ies in the sample, as follows. Photometry and Re values 
are from de Vaucouleurs fits to the light profile. The 
radii and a values are downloaded from the SDSS Cat- 
alog Archive Serveip and the apparent ph otometry from 
the N YU Value-Added Galaxy Catalog (jBlanton et al.l 
120051 ). The photo metry is subsequent ly corrected for 
Galactic extinction (jSchlegel et al.|[l998[ ). K-corrected to 
the F-band at z = usin g the IDL code kcorrect v4.1.4 
(jBlanton &: Roweisll2007() . and converted to absolute lu- 
minosities assuming a standard ACDM cosmology with 
Oa = 0.7, Qm = 0.3, and ho = 0.7. Both L and 
Re have been corrected for known problems with the 
SDSS pipeline sky-subtraction around bright galaxies, 
as described in Paper I. The surface brightness le is 
then computed as /g = Lv/^ttRI. All surface bright- 
nesses and luminosities used in this work are measured 
in the T^-band. The a values are aperture-correct e d to a 
cons tant Re/S apert u re follo wing [jorgensen et al.l (|1995l ) 
and iBernardi et al.l (|2003a[ ). Spectra for the sample 
galaxies were downloaded from the SDSS Data Archive 
ServeiQ. 

The individual SDSS spectra of our galaxies typically 
have S/N r^ 20 A^^, while accurate stellar population 
analysis requires S/N > 100 A" ^ (jCardiel et al.l 11998^. 
It is possible to make estimates of galaxy stellar pop- 
ulation properties f rom individual SDSS spectra (e.g., 
iGallazzi et al.l[2005l ) but we have chosen to take a differ- 
ent approach by defining bins that group together "sim- 
ilar" galaxies, then stacking the spectra of all galaxies in 
each bin to produce high S/N mean spectra that track 
the average properties of galaxies throughout the binning 
space. The ideal S/N > 100 A~^ requires ^ 25 galaxies 

* This emission cut excludes known star-forming galaxies and 
also LINER hosts. The LINER hosts can be considered quiescent 
systems, in the sense that they do not contain ongoing star for- 
mation, but they h ave systematically younger stellar populations 
llUraves et al.|[2"007l ) and have larger corrections to H/3 for emission 
infill. We will return to LINERs in a future paper. 

^ http://cas.sdss.org/dr6/en/ 

^ http://das.sdss.org/DR6-cgi-bin/DAS 
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Figure 1. Sorting and binning galaxies in 3D FP space, (a) An edge-on view of the best fit to the FP, using a least-squares fit of le onto 
Re and a, which gives /{a, Re) = l.lGlogcr — 1.21 log i?e + 0.55. We divide the FP into three slices based on surface brightness residuals 
(Alog/e) to define low-surfacc-brightncss, midplanc, and high-surface-brightness slices, (b) Within each slice, galaxies are further divided 
into 6 bins in log cr and 5 bins in log Re , as illustrated. Diagonal shading indicates the six bins that do not contain enough galaxies to 
produce robust mean spectra and are therefore excluded from our analysis. Vertical shading indicates a bin where stellar population ages 
cannot be reliably measured (see Paper II). 



per bin, depending on the S/N of the individual spectra 
that go into the composite. 

In the analysis presented here, we use two different 
binning strategies to explore mass-to-light ratio varia- 
tions along the FP and through its thickness. The first 
binning scheme is the one used in Paper II, which sorts 
and bins galaxies in the 3D parameter space defined by 
a, Re, and A/eQ This parameter space is ideal for study- 
ing trends in galaxy properties along a (nearly) face-on 
projection of the FP itself. The edge-on projection of 
this plane is shown in Figure [T^. We divide FP space 
into three slices based on residuals from the FP in the le 
dimension (Alog/e). These include a "midplane" slice, 
as well as high- and low-surface-brightness slices above 
and below the plane. We then define a 6 x 5 grid in log a 
and logi?e within each shce (Figure [It). Bin widths in 
logo-, log Re and logA/^ (0.09 dex, 0.20 dex, and 0.10 
dex, respectively) are substantially larger than the un- 
certainties in the measured quantities (typically ^ 0.04 
dex, ~ 0.02 dex, and ~ 0.03 dex, respectively) so that 
the assignment of a galaxy to a given bin is relatively ro- 
bust. We exclude bins that include fewer than 10 galax- 
ies and/or have composite spectra with S/N < 50 A~^ 
(shaded regions in Figure [Tb). 

In §5.2l wc examine trends through the thickness of the 
FP. For this section, we bin galaxies in a 2D binning 
space defined by a and A/e only. In Figure [U we show 
the distribution of the sample galaxies in cr-A/e space, 
which is essentially a cross-section projected through the 
thickness of the FP. We define a 6 x 5 grid in this space 
and use it to sort galaxies into bins through the FP cross- 
section. Bins in a are chosen to be consistent with the 
3D binning scheme of Figure [TJ Bins in A/e span the 
majority of the data; the bin width is chosen such that at 
high cr, even the outlying A/e bins contain > 20 galaxies 
and produce stacked spectra with S/N > 100 A~^. This 
binning scheme means that galaxies with a variety of 



"^ For brevity, we will use the terms cr, Re and Ale throughout 
this paper to refer the the various dimensions of FP space. It should 
be understood that all bins are defined in log space (i.e., logcr, 
log Re , and A log le ) and that we always measure each quantity in 
logarithmic units. 
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Figure 2. Sorting and binning galaxies through the thickness of 
the FP. The quantity A log le is defined as in Figure [T^. We sort 
the sample into the same 6 bins in a as in Figure [TJ], and further 
into 5 bins in A log le ■ 

Re values are combined together in each bin. However, 
as shown in Paper II, the stellar population properties 
and therefore the star formation histories of quiescent 
galaxies do not depend on Re at fixed a. 

With A/e defined in this way, there is a direct cor- 
respondence between le variations and variations in 
Mdyn/L. RecaU that Mdyn/L oc a'^Re/hRl oc a'^/hRe- 
Variations in le are defined at fixed points in a and 
Re- Thus variations in le at fixed values of a and Re 
correspond to variations in Mdyn/L, such that A/g oc 

^{Mdyn/L)-\ 

Once the galaxies have been sorted into bins, the spec- 
tra of all galaxies in each bin arc stacked together using 
an algorithm that masks areas around bright sky lines 
and rejects highly deviant pixels. We measure the full 
set of Lick indices in the resulting mean spectra, then use 
them to model the mean stellar population properties of 
the galaxies in each bin using the single burst popula- 
tion models of ISchiavonI (120071) and t he code EZ_Ages, 
described in [Graves fc SchiavonI (|2008D . Sections 4.2-4.4 
of Paper I give a detailed description of this process. 

There is one significant difference between the data 
presented here and those of Paper II: here the stacked 
galaxy spectra have been corrected for small quantities 
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Figure 3. Median values of cr and M^y^/L for each bin of galax- 
ies (3D FP binning). The listed values of Re are in kpc. The 
solid line shows a least-squares fit of M^y^/L onto a for galaxies 
on the FP midplane, with the slope of the relation indicated in 
the lower right corner. If the FP and the virial plane were equiv- 
alent, M^yn/L would be constant for all galaxies. Instead, there 
is a strong increase in M^y„/L as a increases (the ^^a tilt" of the 
FP) and a substantial spread in M^y^/L at fixed a and Re (the 
"thickness" of the FP). 

of emission infill in the H/3 absorption line used to deter- 
mine ages, as described in Appendix [Xj The quantity of 
emission infill is similar in the various stacked spectra. 
Therefore the main effect of this correction is to system- 
atically lower the derived ages by ^ 0.f2 dex and to raise 
[Fe/H] values by ^ 0.06 dex compared to Paper II. 

3. PHYSICAL EFFECTS THAT CONTRIBUTE TO Moyn/L 
VARIATIONS THROUGHOUT THE FUNDAMENTAL 
PLANE 

As discussed above, both the tilt of the FP and its finite 
thickness represent departures from the simple assump- 
tion of constant Mdyn/ L. This is illustrated in Figure 
[3] for the galaxies in our sample, which span almost an 
order of magnitude in Mdyn/L. The systematic variation 
of Mdyn/L with (T is a 2D projection of the 3D "tilt" of 
the FP. This a-component to the FP tilt has a slope of 
Mdyn/L ex (tO-95 Jj^ |-j^g data presented here, based on a 
least squares fit that weights each binned data point by 
the number of galaxies in the bin. The tilt of the FP 
is typically parameterized as a function of Mdyn, but it 
is sometimes useful to discuss the a component of the 
tilt separately (as in ij5.1[ where we show that the stellar 
population contribution to the tilt is purely a function of 
a). Throughout this paper, we will use the term "cr tilt" 
to discuss mass-to-light variation along the FP as a func- 
tion of a and the term "total tilt" to describe variation 
along the FP as a function of Mdyn- 

In addition to the a tilt shown in Figure [3l Mdyn/L 
also varies significantly at fixed cr. This variation is dom- 
inated by differences in A/g through the thickness of the 
FP (different symbols in Figure[3l), although Mdyn/L also 
varies with R^ at fixed a (different colors in Figure [3]). 

The FP appears to be th e same in all environ- 
ments fKochanek et all 120001: Ide la Rosa et all 120011: 
Ide Carvalho et al., ,2003t iReda etal., ,2005k but see 



iBernardi et al.l [2003b[ ) and inclu des both elliptical and 
SO galaxies on t he same plane (iJorgensen et al.l 119961 : 
iFritz et al.ll2005[) . The slope of the FP evolves with 
redshift, consistent with the pas s ive ag i ng of the stellar 
popu lati ons (e.g.. iKelson et al.l 119971: iJeirgensen et al " 
1999|: IBernardi et al.l l2003bl: Ivan de Ven et al 



20031 iHoTden et al.l 120051: iJoreensen et al.l 120061: 



van der Marel k, van DokkumI I2007D or with pas 



sive evoluti on plus a smal l quan tity of residual star 
formation (jGebhardt et al.l l2003| ). hnt t he thickness 
of the FP is constant out to z ^ 0.6 ([Kelson et all 
[l997Hj0reensen et al.l[l99l iTreu et al.ll200lD . Thus the 
tilt of the FP appears to be time-dependent but the 
thickness of the FP does not, suggesting that these two 
departures from constant Mdyn/L may be driven by 
different mechanisms. 

A number of different physical processes can contribute 
to the variations in Mdyn/L that produce the tilt and 
thickness of the FP. These include differences in the dy- 
namical structure of galaxies, variations in the ratio of 
dark matter to stellar mass, variations in the intial mass 
function (IMF) with which galaxies form their stars, and 
difference in the stellar mass-to-light ratios due to the 
star formation rates, ages and metallicities of the galax- 
ies. 

To elucidate the contributions of each of these physical 
processes, we separate the quantity Mdyn/L into four 
components as follows: 



M, 



dyn 



M, 



dyn 



Mt, 



M, 



L 



Mt, 



M^ M^JMF 



X— ^, (1) 



where Mdyn is the dynamical mass estimate for the 
galax y given by Mdyn = 5a^Re/G fe.g.. lCappellari et al.l 
120061 ) ■ Mtot is the total mass of the galaxy inside Re, M^, 
is the true stellar mass of the galaxy inside Re, Mi,j]\jp 
is the estimated stellar mass inside Re assuming some 
constant IMF, and L is the F-band luminosity inside Re 
(all quantities include L are projected within Re). With 
these definitions, the components of equation [T] represent 
the following physical properties: 

• Mdyn/Mtot is the "dynamical structure term." 
This accounts for any discrepancies between the 
simple dynamical mass estimator and the total pro- 
jected mass within Re. 

• Mtot /Mi, is the "dark matter term." It accounts for 
differences between the total projected mass within 
Re and the projected luminous stellar mass within 

Re- 

• Mi,/Mi,jMF is the "IMF term." It accounts for 
differences between the true projected stellar mass 
inside Re and the projected stellar mass computed 
with an assumed initial mass function. If the true 
IMF matches the IMF used in the stellar popula- 
tion modelling, Af^/Af^jj\fi? = 1. 

• Mi,jMF/ L is the "stellar population term." It 
accounts for the stellar mass-to-light ratio, given 
an assumed IMF. The stellar population term can 
be estimated from the mean stellar age, the mean 
metallicity, and the star formation history of the 
galaxy. 
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Unfortunately, not all of these components can be de- 
termined observationally. For the SDSS data used in 
this analysis, only Mdyn, M^^jme, and L are measurable 
quantities. In practice, M* can only be robustly deter- 
mined for resolved stellar systems. M^, values quoted for 
unresolved stellar systems always rely upon stellar popu- 
lation modelling plus the assumption of a single, constant 
IMF. They are therefore really values of Mi,jMF, not a 
true measurement of M^, . 

Direct measurements of Mtot are only possible where 
an independent mass estimated can be derived. Re- 
cent examples include resolv ed kinematic observatio ns of 
nearby early type galaxies (iCapDellari et al.l 12006) and 
strong galaxy-galaxy lenses (iBolton et al.l 120071 120081 : 
iKoopmans et al.ll2009D which show that the simple dy- 
namical estimator does an excellent job of reproducing 
the galaxy masses derived from kinematics or lensing. 
These studies suggest that M^yn = Mtot, with scat- 
ter of 0.06-0.07 dex standard deviation around the one- 
to-on e relation (cf. iCappellari et al.l [20061 IBolton et al.l 
|2007() . This level of variation in Mdyn/Mtot is clearly 
inadequate to explain the variation in Mdyn/ L in Fig- 
ure [3] These observational results are further sup- 
ported by numerical simulations which show that galax- 
ies with orbital anisotropics large eno ugh to move them 
off the FP are dy nami cally unstable (iCiqtti et al.l 119961 : 
iCiotti fc Lanzonilll997l: iNipoti et alll2002D . 

Assuming that M^yn = Mtot and bearing in mind that 
Mi, is not a measurable quantity in unresolved stellar 
systems, equation [1] simplifies to 
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where both Mdyn/Mi^jj^jF and M^^jmf/L are measur- 
able quantities. Variations in M^^jmf/L are due to 
known stellar population effects, while variations in the 
quantity Mdyn/Mi,jMF can be produced by variations 
either in the IMF or in the dark matter fraction inside 
i?e- Hereafter, we will refer to the Mdyn/Mi,jMF term 
as the "dark matter/IMF" contribution with the under- 
standing that variations in this term can be attributed 
to either effect, or to both in combination. 

Many previous studies have indicated that the 
stellar population term {Mi,jMF/L) contributes 
not more than 1/2 of the ob serverd F P tilt (e.g., 
Pahre et al.l 119951: iPrugniel fc Sim ien 19 9'7t iPahre et alJ 



1998t IScodeggioetaLJ 119981: iMobas her et al.l 
Padmanabhan et al.l [ 2004 : 



iTruiillo et al 



1999 : 



2004' 



Gallazzi et al.l 120051: iJun fc Iml [2008t iLa Barbera et al 



2008HHvde fc Bernardill2009l . but see also lAllanson et al 



20091 who find that using a single burst star formation 
history produces enough variation in Mi,jMF / L to 
explain the full FP tilt). This then requires vari- 
ations in M(iyn/Mi,jMF d ue to variations in e ither 
the IMF (a. s prop osed by iRenzini fc Ciottil 119931 and 
iChiosi et al.l Il998f )_ or in the dark matter fraction 
inside Rp 



iM^ 



Renzini fc Ciotti '1993': ' Pahre et al.l 
200 It iPadmanabhan et a l. 2004; :Gallazzi et all I200E 



Hvde fc Bernardil I2009f ) to produce the observed tilt of 



the FP. 

The case for dark matter variation has recently gotten 
a boost from simulations of galaxy mergers that include 
gaseous dissipation. Several authors have shown that 



the degree of dissipation during mergers determines the 
dark matter fraction within R^. in merger remnants (e.g., 
iRobertson eral|[200l [Pekel fc Cox|[2006[ ). Under realis- 
tic conditions in which mergers of more massive galaxies 
have decreasing amounts of gas, such si mulatio ns can re- 
produce the observed t ilt of the FP (iKobayashi 2005t 
Robertson et al.l [20p [Dekel fc Coxl 120061: iCoviilgto^l 



2008HHopkins et al.ll2008aD . A limited quantity of fur- 



ther dissipationl ess mer ging should pr eserve the tilted FP 
jCapelato et al." '1995": Dantas et al.' '2003': 'Nipoti et aU 
[2003; Boylan-Kolchin et al. 2005; Robertson et al. 200^ 
It seems therefore that, while stellar population effects 
may contribute to the tilt of the FP, further variation 
in the dark matter/IMF term are required. This paper 
shows that these two contibutions to the FP tilt rotate 
the plane around different axes in the 3D space, such 
that M^jMF oc O-0-30, while Mdyn/M^jMF oc M°yf. We 
also show that dark matter/IMF variations are required 
to produce the thickness of the FP: stellar population 
effects contribute some, but the dominant component of 
variation must come from differences in the dark mat- 
ter/IMF term. 

4. MEASURING Mi,jMF/L 

As discussed in section [31 the quantities Mdyn, 
Mi^jMF, and L are the only components of equation 
[1] that can be determined for the SDSS data. Of 
these. Mdyn and L are relatively direct measurements 
while Mi,jMF, the stellar mass estimate, comes from 
L X Mi^jj^ip/ L, where M^^jmp/L is determined through 
stellar population modelling. In Paper II, we measured 
the mean age, metallicity, and abundance patterns for 
the stacked spectra in 3D FP space. With these quan- 
tities determined, there are three remaining factors that 
affect Mi,jMF/L. These are: 

1 . The stellar population model code used to compute 

Mi,jMF/L. 

2. The library of star formation histories used in the 
models. 

3. Whether the stellar population modelling is done 
by fitting individual spectral absorption lines, or 
by fitting the broad-band galaxy spectral energy 
distributions (SEDs). 

This section addresses how each of these factors affect 
the determination of AI^jmf/ L, in order to motivate our 
choice of models and to assess the extent to which our 
final conclusions depend upon this choice. Qualitatively, 
our conclusions turn out to be independent of the exact 
choice of stellar population model and star formation his- 
tory used to determine M-^jmf/L, although quantifying 
t hese relationships does depen d on the choice of model. 
iLonghetti fc Saraccol (|2009[ ) have recently addressed 
the first issue listed above using a suite of stellar pop- 
ulation models to predict Mi,jMF/L based on opti- 
cal and infra-red colors. They examine the depen- 
dence of Mi,jMF/ L on age, metallicity, the IMF, and 
the choice of stellar population model codes. They 
ex plore several popu l ar mo del c odes inclu d ing t hose 



ex plore several popu l ar mo del c odes inclu d ing t nose 
of iBruzual fc CharlotI ([200l and iMarastonI ([2005[) . as 



well as P EGASE (iFioc fc Rocca-VolmerangeJll997D ! and 
GRASIL (|Silva et al.l ll998D) and find that the different 
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Figure 4. Age and [Fe/H] values measured in tiic stacked spectra (3D FP binning). Multiple symbols with the same size and color 
represent different values of Re; only bins with small errors (A[Fc/H] < 0.05 dex and AAge < 15%) are shown. Galaxies with higher values 
of (T typically have older ages and higher [Fe/H] than galaxies with low values of a. At fixed a, galaxies with high (low) surface brightness 
have younger (older) ages and higher (lower) [Fe/H] than galaxies at the same a that lie on the FP midplane (c.f., Paper H). Lines of 
constant M^,i]^jp/L (a) and lines of constant g — r color (b) from BC03 models are overplotted as the gray lines. From a given value of 
age and [Fe/H], the corresponding value of Mi,jMF/L and g — r can be estimated. Mi,jMF/ L depends strongly on age but only weakly 
on [Fe/H] , while g — r depends strongly on both age and [Fe/H] . The anti-correlated variations in age and [Fe/H] at fixed cr move galaxies 
nearly along lines of constant g — r but spread the galaxies over a range of Mi,jMF/L. 



codes produce relatively consistent M^^jnip/L predic- 
tions. 

If the ages and metalliciti es of the galaxies are known, 
iLonghetti fc Saraccol (I2009D find that the largest varia- 
tions in Mi,jMF/L determinations are due to the choice 
of IMF, w ith a zero point shift of ~ 0.26 dex between the 
ISaloeteH ([Ml) and lChabrieJ (12001 IMFs and of ~ 0.21 
dex between the Salpeter and iKroupal ()2001[ ) IMFs, in 
the sense that the Salpeter IMF predicts larger values 
of Mi,jMF/L than either Chabrier or Kroupa. Relative 
values of M-i,j]\[f/L between different galaxies are ro- 
bust to this effect, as long as the IMF is the same for all 
galaxies. 

In the following sections, we examine the effects on 
Mi,jMF/L of various star formation histories, and com- 
pare SED-based estimates of Mi,jMF I L with those de- 
rived from spectral absorption features. All of the mod- 
els presented here are based on isochrones computed for 
scaled-solar element abundance patternsO 

® Since the galaxies presented here typically have super-solar 
[Mg/Fe], these are not an ideal match to the data. Mg-enhanced 
isochrones are shifted to cooler temperatures than scaled-solar 
isochrones, but a complicated interplay between different abun- 
dance effects makes it challenging to predict ho w various non-solar 
abundance patterns will affect M^- imf/ L (e.g.. lDotter et al1l2007l : 
ILee et al.l2010l ) . Understanding these effects in detail may well rep- 
resent the next major step forward in stellar population modelling. 



4.1. Single Burst Models for AI-^j]\[f/L 

As a reference point, we begin by computing our own 
values of M^^jmf/L using simple single-burst stellar pop- 
ulation models (SSPs). Having mapped out the varia- 
tions in mean age and metallicity through 3D FP space 
in Paper II, we are now in a position to study how 
these variations contribute to the observed variations in 
Mdyn/L. The modelling process is illustrated in Fig- 
ure m which shows the values of mean, light-weighted 
age versus [Fe/H] determined from the stacked galaxy 
spectra, as described in Paper II. We have overlaid 
lines of constant stellar population mass-to-light ratio 
{Mi,jMF/L, panel a) and constant g — r color (panel 
b) de rived from the single burst stellar population mod- 
els of iBruzual fc CharlotI ()2003L hereafter BC03) with a 
Chabrier IMF. These curves make it possible to read 
off model values of Mi,jMF/L and 5 — r for each set 
of stacked spectra considered here0 

The different colors and symbols encode the a values 
for each stacked spectrum, while the symbol size encodes 

^ The BC03 models are computed for solar abundance ratios 
only, which is not an accurate reflection of the abundance pattern 
of our sample galaxies. Nevertheless, we show in Appendix IB] that 
the BC03 models, when chosen to match the SSP ages and [Fe/H] 
values from our stellar population modelling, are able to accurately 
reproduce the observed galaxy colors. 



Mass-to-Light Variations in 3D FP Space 



Ale, as indicated. The plot includes the various bins 
in Re but these are not differentiated in the figure, as 
neither age nor [Fe/H] depend on Re (see Paper II). The 
figure illustrates the main conclusion of Paper II: that 
the stellar populations of quiescent galaxies form a two- 
parameter family. The two dimensions of this family are 
such that galaxies with higher values of a tend to have 
older ages and higher [Fe/H] than galaxies with low a, 
while at fixed a galaxies with lower surface brightnesses 
have older ages and lower [Fe/H] than those with higher 
surface brightnesses. 

It is interesting to note that, although galaxies at fixed 
a span a significant range in age and [Fc/H], this varia- 
tion moves them almost exactly along lines of constant 
g — r color, as shown in Figure |4)d. Thus, a process that 
attempts to use galaxy colors to model stellar popula- 
tions will not successfully distinguish between the differ- 
ent populations at fixed a. This is due to a nearly perfect 
alignment of the age-metallicity degeneracy in color with 
the measured anti-correlated of age and [Fe/H] in galax- 
ies at fixed (T. However, lines of constant M^^jme/L in 
Figure 13^ follow shallower trajectories compared to lines 
of constant g — r color, such that galaxies at fixed a span 
a range of M^^jMr/L values but not a range of colors. 
Using the full 5-band SDSS photometry does not give 
significantly improved Mi,j]\jp/ L values, as all optical 
colors follow similar trajectories in age-[Fc/H] space F^ 

The models shown here are based on single burst stel- 
lar population models (SSPs), which are almost certainly 
not accurate descriptions of the star formation histories 
of galaxies. However, such simple models are more trans- 
parent to interpret, and they provide a good reference 
point for considering models with more complicated star 
formation histories. Appendix [B] shows that the g ~ r 
colors derived from Figure |3Jd do an excellent job of re- 
producing the observed galaxy colors, which gives us con- 
fidence that the derived values of Mi^jMp/L are likewise 
reasonable. 

4.2. Models for Ali^jMp/L with Complex Star 
Formation Histories 

For galaxies with old stellar populations, the differ- 
ences between single burst values for Mi,jMF/L and 
those from extended star formation histories are small. 
However, galaxies with younger mean light-weighted ages 
may harbour a bright young sub-population obscuring a 
large population of older stars that contribute substan- 
tially to Mi, but very little to the integrated galaxy light. 
For these galaxies, the single burst estimates for younger 
galaxies are almost certainly too low. 

This is illustrated in Figure[S^, which plots AI-^jmf I L 
against the H/3 absorption line strength^ The figure 
compares single burst models (dashed line) to models 
with continuous star formation (solid line) that begins 
at z = 10, proceeds at a constant star formation rate, 
then is abruptly truncated at various timesQ These 
star formation histories are indistinguishable for values 

^^ There is a slight steepening of the Hnes of constant color as 
colors become bluer (e.g., lines of constant u — g are steeper than 
g — r) but even the reddest i — z color is not nearly as flat as the 
lines of constant Mi,j},ip/L in Figure|4] 

^^ H/9 is the major age discriminant used in the stellar popula- 
tion modeling of Paper II. 

^^ In order to isolate the effect of age distribution in the stel- 



of \og{Mi,jMF/L) > 0.35, corresponding to single burst 
ages > 7 Gyr or constant star formation rates with 
truncation times it > 3 Gyr ago. At lower values of 
AI*,imf/L, single burst models can substantially under- 
estimate Mi,jMF/L if the true star formation history is 
more extended. 

In Figure EJ;, we show H/3 versus M^^^mf/L for a set 
of exponentially declining star formation rates with e- 
folding time r (so-called "r models"). For each r, we 
compute a set of models in which the onset of star forma- 
tion occurs at different times, spanning a range formation 
redshifts with 0.25 < Zf < 40. However, the majority of 
the models with r > 0.5 imply z = Q star formation rates 
that are too large for our quiescent galaxy sample. These 
are excluded from the set of viable models, as explained 
below. 

The sample presented here includes only galaxies with 
Ha emission-line fluxes below a 2ct detection threshold. 
We take the 2(T error value as an upper limit on the intrin- 
sic Ha flux for each galaxy, the n convert this lim it into 
a star formation r ate following iKennicuttl ()1998l equa- 
tion 2). Using the IGallazzi et al.l (j2005[) rneasurements 
of A'hjMF for each galaxy (explored in greater detail in 
> j4.3|) . we compute an upper limit to the specific star for- 
mation rate (SSFR) for each sample galaxy. The vast 
majority (96%) of our sample galaxies have upper lim- 
its on their SSFRs that are below 0.1 Mq yr-7lO"M0. 
We therefore use this value to restrict the set of r models. 

In Figure [5]:, we show only models that have SSFRs 
at z = that are less than 0.1 Mq yr"VlO"^0- Low 
values of r lead to models with short e-folding times, 
which closely resemble single burst star formation his- 
tories. At high values of r, only models in which star 
formation begins at very early times have low enough 
SSFRs today to be included in the sample. The net ef- 
fect is that the Mi,jMF/L values for a given observed 
H/3 line strength are nearly identical to the single burst 
value, regardless of the values of r or zj. For t- model 
star formation histories, single burst models in fact give 
a good approximation of the true value of M-^jmf/L. 

However, models which include an abrupt truncation 
of star formation can have substantially different values 
of Mi,jMF/L for a given H/3 line strength. In Figure 
[SJi, we modify our library of t models to include models 
that follow exponentially declining star formation histo- 
ries that are suddenly truncated 1 Gyr before the time 
of observation. These models show some variation in H/3 
line strengths for a fixed value of M-i,jmf/L, depending 
on the star formation history. The spread is modest for 
galaxies with Mi,jMF/L > 0.35 but is larger for galaxies 
with lower M-^jmf / L. 

Based on this suite of results, we infer that M-^jmf I L 
values above 0.35 are relatively insensitive to the details 
of the galaxy star formation histories and are therefore 
rehable. In contrast, values of M^^jmfIL below 0.35 
may be off by 0.1-0.15 dex and are more likely to be 
skewed too low for the youngest galaxies with the lowest 
M^jmfIL^ 



lar populations, the constant star formation rate models, as well 
as the more complicated r-models considered in this section, are 
computed at Solar mctallicity. 

^'^ It should be noted that these statements apply to smoothly 
varying stax formation histories. Non-continuous (i.e., bursty) star 
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Figure 5. A comparison of Mi,jMp/L values from models witii different star formation histories. All models are BC03 models with a 
Chabrier IMF and solar abundances, (a) Simple star formation histories: the dashed line shows M^^jmf/L and H/3 values from single 
burst star formation histories, while the solid line shows values from models with constant star formation that begins at ^ = 10 and 
truncates abruptly at various times. Differences in M^jj^jp/L due to star formation history become relevant for \o^{M^,iMp/ L) < 0.35, 
which corresponds to iq„enc/i = 3 Gyr ago. (b) Models with a frosting of metal-poor (Z = 0.0001) stars, with fractions as indicated. 
Small, very metal-poor populations produce slightly higher M^jmf / L for a given value of H/3. (c) r-models with 0.1 < t < 4.0 Gyr 
and a range of initial formation redshifts (0.25 < Zf < 40). Only models with SSFRs consistent with our sample galaxies are shown 
(< O.IMq yr~^/lO^^M0). These all have mass-to-light ratios consistent with the single burst models, (d) r-models that have been 
abruptly truncated 1 Gyr before the time of observation. This family of models produces a limited range of H/9 values for a given value of 
M*,imf/L, particularly for log(Mj,_/MF/i) < 0.35. 



Finally, IMaraston et al.l (|2009f ) have suggested that 
early type galaxies may contain a small sub-population 
of very metal-poor stars. They find that models for 
single-metallicity old stellar populations cannot match 
the color evolution of luminous red galaxies in the SDSS 
from z = 0.7 to z — 0.1. They demonstrate that adding 
a frosting (3% by mass) population of low metallicity 
{Z = 0.0001, i.e., [Fe/H] = -2.2) stars makes it pos- 
sible to match the observed color evolution. In Fig- 
ure [SId, we illustrate the effect of such a population on 
M*.imf/ L and H/?. The 3% fr osting population sug- 
gested by IMaraston et al.l (j2009() results in M-i,jmf/L 
values that are < 0.05 dex higher for a given observed 
values of H/3, which is comparable to the variation pro- 
duced by the r-model star formation histories considered 
here. 

formation histories can also produce substantially different values 
of Mi,jMF/ L. However, in this analysis, we are examining stacks 
of many dozens of spectra, which means that the assumption of 
continuous star formation histories is reasonable. This is one ben- 
efit of stacking. 



This leads to a key point. In all cases, the single hurst 
models give a lower limit on M-^jmf/L. Single burst 
models in fact give us a powerful tool because we can un- 
derstand the biases in measurements based upon them. 
This gives us a benchmark against which to compare 
other, more complicated determinations of Mi,jMF/L for 
which the biases are less clear. 

4.3. Mi,j]\jp/L from. Spectral Absorption Lines vs. 
Optical Colors 

Figure [5] compares M-i,jmf/L from four different 
sets of authors, two of which are co mputed based 
on fi ts to spectral absorption feat ures (jGallazzi et al.l 
120051 and iKauffmann et all l2003bf ) and two of which 
are based on optical colors or multi-band photome- 
try (jBlanton fc RoweisI 120071 and iBell et al.l l2003[ ). to 
our benchmark single burst models. All four of these 
Mi,ji\ip/ L measurements are available for the individual 
SDSS galaxies used in our sample. Each papers' version 
of Mi,jMF /L uses slightly different stellar population 
models and forms for the IMF. Wc have corrected the zc- 
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Figure 6. (a-e) Various estimates of M^, imfI L for the galaxies in our sample as a function of a. Galaxy bins on the FP midplane are 
shown in black, while those on the low- and high-SB slices arc shown in gray. Panels ( a -c) sho w values derived from spectral absorption 
features, including measurements from (a) 'Gallazzi ct al. (2005), ih) IKauffmann et al.l Il2003b|). and (c) our single b urst valu es. Panels 
(d-e) show values based on fits to multi-band photometry, following the prescriptions of (d) iBlanton &: Rowe3 II2007I ) and (e) IBell et all 
1120031 ) . In each panel, solid lines show least squares fits of Af* /j\/i?/L onto a for the FP midplane galaxies only, weighted by the number of 
galaxies in each bin. The slopes are indicated in each panel. The dotted line at M^ imf/L = 0.35 marks the point below which differences 
in the star formation history significantly affect the determination of M*^/jvfi?/L (see section [4.211 . (f) The values of M^y^^/L from Figure 
[3] are reproduced for comparison. Regardless of the method used to measure Mi^jMp/L, it clearly cannot provide either the full slope of 
^dyn/L with (T (the "(T tilt" of the FP) nor the spread in M^y^/L at fixed a (the "thickness" of the FP). 

eled in a comparable way using a similar library of star 
formation histories. Gallazzi et al. use a larger number 
of spectral absorption features in their fits; in particular, 
they include a number of redder features at ^ 5000 A. 

Figure IH] shows that the two spectroscopic stellar 
mass estimates agree very w ell for high-cr galaxies, 
with IKauffmann et al] ()2003b( ) finding somewhat lower 
Mi^jMp/L for low-cr galaxies. This may be because their 
measurements only use absorption lines around 4000 A 
and are therefore more strongly influenced by young sub- 
populations of stars than arc the Gallazzi et al. estimates, 
resulting in underestimates of M^^^mf / L. This differ- 
ence produces a stronger slope to the \og[Mi,jMF/L)- 
log (7 relation (gray line) in Kauffmann et al. than in 
Gallazzi ct al. 

These values are compared to our single burst values in 
panel c, which are similar to the Kauffmann et al. values 
for most galaxy bins. The largest discrepancies are for 
galaxies with low a and low M^ji^ip / L^ where the single 
burst estimates arc significantly lower than Kauffmann 
et al. These are precisely the galaxies for which we expect 



ropoin ts of the IKauffmann et all (j2003b[ ) and IBell et ahl 
(|20C)3|) values to match t he Chabrier IMF assum ed by 
the lGallazzi et al.l (|2005D . IBlanton fc RoweisI (|2007[ ). and 
our sin gle burst models. 

The (Kauffmann et al.l (|2003b( ) and iGallazzi et aD 
(j2005| ) stellar mass estimates for the individual SDSS 
galaxies were downloaded from the der ived data catalogs 
from S PSS studies at MPA/JHlQ. IKauffmann eTaLl 
(|2003bD determine M^,jmf/L b y fitting the Dn4000 
jBalogh et al.l 119990 and H^^ (jWorthev fc Ottavianil 
Il997t ) spectral features to a library of models. Their 
model star formation histories comprise a set of ex- 
ponentially declining star formation rates with varying 
T, a range of total metallicitics (always with the solar 
abundance pattern), and superposed random bursts of 
star formation at late times. We have corrected the 
Kauffmann et al. stellar mass measurements from SDSS 
Petrosian magnitudes to de Vaucouleurs magnitudes to 
match the other ma ss estimates presented here. The 
IGallazzi et al.l (|2005D stellar masses (panel a) are mod- 

^"^ http://www.mpa-garching.mpg.de/SDSS/ 
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the single burst approximation to break down, resulting 
in single burst M^, imf / L values that are too low (section 
[4:2]) . 

The photometry-based determinations of M-^jmf I L 
look different from all the spectroscopic values, as shown 
in Figures IBJi and[6t. The stellar masses shown in Figure 
[5Ji are computed using M. Blan ton's IDL code kcorrect v. 
4.1.4 (jBlanton fc RoweisI [20071 ) . This code matches the 
SDSS ugriz photometry (K-corrected to z = 0) to a lin- 
ear combination of a set of basis templates that have 
been determined to span the space of SDSS galaxies. 
The stellar masses in Figure [BJ; are computed using the 
B — V galaxy colors (determined using kcorrect) and the 
conversions be tween galaxy col or and Mi, imf I L given 
in Table 7 of iBell et all pOOSD . The BeU et al. con- 
versions are based on fits of the full ugrizK combined 
SDSS and 2M ASS photometr y for SDSS Early Data Re- 
lease ((Stoughton et al.l I2002D galaxies to a set of mod- 
els spanning a range of star formation histories. In both 
cases, the colors are computed from photometry that has 
been matched to the SDSS 3" spectral fiber aperture in 
order to provide consistency with the spetroscopically- 
derived values of Mi^jMF/L (see appendix IBJ) . These 
photometrically-derived values of M-^j]\[f/L show much 
weaker variation with a and much smaller spread at fixed 
a than the spectroscopically-derived values. 

The increased spread in the spectroscopic values can 
be understood from Figure HId. As discussed in section 
14.11 the lines of constant 5 — r in the models of Fig- 
ure HJd run very nearly parallel to the lines of constant 
(T in the data. The spectroscopic data are thus capable 
of differentiating between different star formation histo- 
ries and different Mi^jMp/L values at fixed a, which the 
photometry-based values cannot do. 

Throughout the re st of this analysis, we will use the 
IGallazzi et al.l (|2005f ) values for M^^jmf/L because they 
capture the variation in M^^jmf/L at fixed a (unlike 
the photometry-based measurements) but are likely to 
be less biased toward younger sub-populations than the 
[Kauffman n et al.. (,2003bl ) and especially the single burst 
values. 

4.4. Mdyn/L versus A^jmf/L 

As a final point of comparison. Figure |Bf shows Mdyn/L 
as a function of a. It is clear that none of the Mi,jMF/L 
measurements (panels a-e) resemble Mdyn/L. No matter 
which method is used to measure M-i,jmf/L, stellar pop- 
ulation effects cannot reproduce the observed variation in 

Mdyn/L. 

In the first place, none of the M^^jhif/L mea- 
surements have a slope with a that is as steep 
as the slope of the observed Mdyn/L-a relation. 
This means that stellar population effects are not 
enough to explain the tilt of the FP, as has been 
demonstrated by numerou s previo us authors (e.g., 
Pahre et al.|[T99a IPrugnie]_fc Simiei lfTOOl IPahre et alj 



19981 iScodeggio et all [I99& .Mobasher et al.l 1199 



Padmanabhaii et al.ll2004i: iTruiillo et al.ll2004HjInr fc Ir 
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^^ In contrast to our results and other previous studios, 
lAllanson et al.l II2009I) find that the assumption of single burst star 
formation histories predicts a strong enough scaling of M^^jme/L 
with a to account for the entire a tilt of the FP. While our re- 



in the second place, it is clear from Figure [S] that 
M^^jmf/ L variations cannot provide enough spread at 
fixed a to match the observed variation in Mdyn/L. This 
new result indicates that stellar populations contribute lit- 
tle to the thickness of the FP, despite the fact that resid- 
uals from the FP app ear to be correlate d with stellar 
population properties (iForbes et all Il998t IWuvts et al.l 
l2004HGargiulo et al.ll2009l Paper II) This is because the 
stellar population differences through the thickness of 
the FP show an anti-correlation between age and [Fe/H] . 
This anti-correlation, while not exactly parallel to lines 
of constant Mi,jMF/L (see Figure |4]), is in a similar di- 
rection and therefore limits the differences in M^^jmf/L 
that result from the stellar population differences at fixed 

G. 

The marked difference of Figure [6f from Figure [6^-6 
requires that there be variations in Mdyn/Mi, (the dark 
matter fraction) and/or in M^^/Mi^^MF (the IMF) among 
quiescent galaxies. At least one of these quantities must 
vary along the FP to reproduce the observed a tilt of the 
plane. One or both of them must also vary through the 
thickness of the FP. Furthermore, these variations likely 
dominate over stellar population variations both along 
and through the FP. 

5. MASS-TO-LIGHT RATIOS IN 3D FP SPACE 

Having quantified the variation of M-^jmf/L for our 
sample galaxies, as well as the systematic uncertainties 
in Ali,jf,[F / L measurements, we can now examine how 
the stellar population term contributes to the total ob- 
served variation of Mdyn/L throughout 3D FP space. In 
essence, we will be dividing out the variations due to stel- 
lar populations and thereby mapping how the remaining 
term {Mdyn/Mi,jMF) varies over the plane and through 
it. 

Section FS.lff ocuses on the FP midplane, examining how 
Mi,jMF/L and Mdyn/M-i,jMF contribute separately to 
the tilt of the FP. We show that, not only does the tilt of 
the FP require contributions from both Mi,jMF/L and 
Mdyn/Mi,jMF, but that these two components of the tilt 
rotate the plane around different axes in the 3D param- 
eter space. Section 15.21 then maps both these quantities 
in a cross-section through the FP, examining how each 
contributes to the thickness of the FP. We show that 
the thickness of the FP is dominated by variations in 
Mdyn/M-i,jMF, with M-i,jmf/L cont ributing only ~ 22% 
of the variation. Finally, section [531 distills the quantita- 
tive results of this analysis by presenting fitting functions 
for the total mass-to-light variation (Mdyn/Ly), the stel- 
lar population contribution [M-^jmf/L), the dark mat- 
ter/IMF contribution {Mdyn/Mi,jMF), and the associ- 
ated surface mass density of stars (E^,./A/i?) as functions 
of the three FP variables: a-, Re, and A/g. 

5.1. The Tilt of the FP in 3D 

suits agree that single burst models produce substantially stronger 
slopes than those computed with more complex star formation his- 
tories, we find that even single burst models cannot produce the 
entire a tilt of the plane. A comparison between our Figure |6] 
and their Figure 12 suggests that this difference is due to differ- 
ences in the slope of the observed M^y„/L-a relation (they find 
Mayn/L oc O-0-79 while we find May„/L oc o-0.95i0.04-) Qur single 
burst stellar population modelling results agree with theirs; they 
find M-t. IMF /L oc O-0-66 £q]. ^ single burst model, while we find 
M,jMF/i OCaO-67±0.04. 
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Figure 7. The FP tilt contribution from stellar population effects, mapped throughout 3D FP space. Data points in each panel represent 
the sample galaxies divided into three slices in surface brightness, as in Figure[T] The overlaid color contours show variations in M^^jnip/L in 
this space, measured using the values of lGallazzi et aO 1)20051 ) . On the FP midplane (center panel), lines of constant M^-^/mf/^ follow lines 
of constant cr (dashed line); M^ imf/ L is independent of R^. This indicates that M^. imf/L is a function of only a {M^. imf/L = /(o")) 
on the FP. The low- and high-surface-brightness slices show the same effect, with lines of constant M^,j]\fp/L running vertically. However, 
they show zeropoint shifts from the midplane, with the low- and high-surface-brightness galaxies have systematically higher and lower 
J^-k,lMF/ L, respectively. These trends in Mi,jMp/L are as predicted from the stellar population maps in Paper II. 

Figure [7] maps how the stellar population M-^jmf/ L 
component varies over the FP using the bins presented 
in Figure [TJ The three panels show the three slices 
in surface brightness, including the midplane and the 
slices above and below the plane. The overlaid con- 
tours indicate the media n Mi,jMF I L v a lues fo r each bin 
in FP space, using the iGallazzi et al.l ()2005h measure- 
ments. The contours are constructed by plotting the 
median value of Mi^jMF / L for each bin at the point cor- 
responding to the median values of a and Re for the bin. 
The values of M^^jme/L are then linearly interpolated 
between the 6x5 grid of bins in each surface brightness 
slice to produce a continuous map of M^^jme/L along 
each slice. 

Looking first at the FP midplane (center panel), it is 
clear that M^^jMr/L depends on a but is independent of 
Re- Lines of constant M-^jmf I L follow lines of constant 
cr (dashed line). Therefore, we find that M^^jme/L is a 
function only of a on the FP. This is expected from Paper 
II, where we showed that all stellar population proper- 
ties (age, [Fe/H], [Mg/H], and [Mg/Fe]) follow lines of 
constant a on the FP and are essentially independent of 
Re- It therefore follows that the known stellar popula- 
tion contribution to Mdyn/L should also follow lines of 
constant a. 

The low- and high-surface-brightness panels above and 
below the midplane show the same behavior: M-^jmf I L 
again follows lines of constant a and is independent of R^. 
However, the zeropoint of the relation shifts such that 
galaxies in the low-surface-brightness slice have higher 
Mi^jMp/L for a given a than their counterparts on the 
FP midplane. Similarly, galaxies in the high-surface- 
brightness slice have lower Mi,j]\jF/L for a given a. 
This demonstrates that there are Mi^jMp/L variations 
through the thickness of the FP as well, although they are 
relatively weak. We will study these variations through 



the thickness of the FP in more detail in 

With Mi^jMp/L now mapped, we can divide it out 
to isolate the map of Mdyn/Mi,jMF alone. We do this 
by calculating the ratio of the two known quantitites 
Mdyn/L and M-^jmf/L. The resultant map, which 
shows the dark matter/IMF contribution Mdyn/Mi,jMF, 
is shown in Figure [SPi The data and contours are con- 
structed in the same way as Figure [71 

A major result of this paper is illustrated in Fig- 
ure [8] lines of constant Mdyn/Mi^jMF on the FP mid- 
plane (center panel) follow a noticeably different slope 
than the lines of constant Mi,jMF/L in Figure [71 The 
new lines run approximately parallel to lines of con- 
stant Mdyn (dashed line, with Mdyn ex a^Re), indicating 
that on the FP midplane Mdyn/Mi,jMF is a function of 
Mdyn- Similar behavior is seen in the low- an d high - 
surf ace-brightness pan e ls. B oth IGallazzi et al.l ()2005[ ) 
and iHvde fc Bernardil ()2009l ) showed a correlation be- 
tween Mi,jMF and Mdyn that was strong, but this is the 
first time that it has been shown to be the best corre- 
lation, by mapping Mdyn/M^^jMF explicitly in 3D FP 
space. 

By examing the two contributions to the FP tilt in 3D 
space, we have revealed that the two suspected causes of 

1^ The values of May„/M^jMp run from -0.13 to -1-0.50. Neg- 
ative values of M^yn/Mi,jMF are clearly nonsensical if M^y^ = 
^a^R^/G gives a true representation of the total mass and the 
Chabrier IMF is a good approximation to the true IMF in these 
galaxies; a galaxy cannot contain less total mass than its stel- 
lar mass. However, mode st changes to the assumed IMF (e.g., 
ILonehetti &: Saraccol 120091 ) or to the constant used to compute 



M, 



dyn 



will produce zeropoint shifts of the needed magnitude (~ 0.1 
dex). In this work, we have focused on trends through FP space, 
for which the zeropoint is not relevant, rather than on absolute 
measurements of M^,i]^fp/L. It is worth noting however that, as 
computed here, the single burst values of Mi,jMF/L are the only 
ones for which all galaxy bins give Mi,jMF < ^dyn- 
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Figure 8. The FP tilt contribution from varying dark matter/IMF effects, mapped tfirougfiout 3D FP space. Data points in each panel 
represent the sample galaxies, divided into three shces in surface brightness, as in Figure [T] The overlaid color contours show variations in 
^dyn/^-k,lMF in this space. On the FP midplane (center panel), lines of constant M^y^/M^^ji^ip approximately follow lines of constant 
(dashed line) such that M^y„/M^jMF = fi^dyn) on the FP. The two components of the FP tilt therefore rotate the plane 



Mdyn 

around different axes. The low-surface-brightness (high-surface-brightness) slice shows the same effect, but with a zeropoint shift such that 
^dyn/^*,lMF values are higher (lower) than the midplane values. 



FP tilt rotate the plane around different axes. The stel- 
lar population contribution has Mi^j^p/L ~ f{a), while 
the dark matter/IMF contribution has Mdyn/Mi^jj^ip ~ 
g{Mdyn) = g{o'^Re)- This indicates that the tilt of the 
FP cannot be parameterized as a simple power-law func- 
tion of Mdyn or A/* or L or a. Instead, it takes the 
functional form f{a) x g{Mdyn)- 

Figure [5] also illustrates a second major result of this 
paper. Lines of constant Mdyn/M^ji^ip run parallel to 
lines of constant Mdyn in all three panels, but there is 
a zeropoint shift between the slices. The low-surface- 
brightness and high-surface-brightness slices have respec- 
tively higher and lower values of Mdyn/M^,jMF for a 
given a than do their counterparts on the FP midplane. 
The sign is the same as the shift of M^^jMr/L seen in 
the three panels of Figure [T] Both 3D maps reveal that 
the stellar population effects and the dark matter frac- 
tion/IMF effects which contribute the the FP tilt con- 
tribute also to the thickness of the FP as well. Section 
15.21 explores this second trend in more detail. 

Focusing for the time being on the FP midplane alone, 
we quantify the contributions of the two tilt components 
in Figure [9l The left panels show the stellar population 
(upper left) and dark matter/IMF (lower left) contri- 
butions to the a tilt of the FP. Only the galaxy bins 
on the FP midplane slice are included. M-^jmf/L is a 
very tight function of cr. To quantify the stellar popu- 
lation contribution to the a tilt of the FP, we calculate 
the ratio Alog(M^,_7M_F/i)/AlogCT. This is the first of 
several such ratios we will discuss; hereafter we refer to 
them as "A-A" relations. Note that the numeric value 
of Alog(M^,jM_F/i)/Alogcr is equivalent to the quan- 
tity /i in the relation Mi,jMF / L ex a^ . A linear least- 
squares fit, weighted by the number of galaxies in each 
bin, gives Mi,jMF/L ex cr°-^° and scatter of only 0.01 
dex. Meanwhile, Mdyn/M-^jMF is a less tight function 
of CT, with Mdyn/M^jMF oc cr°-^^ and scatter of 0.05 dex. 



Comparing the two slopes shows that dark matter/IMF 
variations contribute twice as much to the a tilt of the 
FP as do stellar population variations. 

We can instead choose to parameterize the FP tilt as 
a function of Mdyn (i-c, the "total tilt" of the FP). In 
this case, the A-A relation between Mdyn/Mi^jMF and 
Mdyn (lower right panel) is significantly tighter than that 
between Mdyn/Mi,jMF and a (lower left panel), with 
Mdyn/Mi,jMF oc M^^^f and scatter of only 0.02 dex. 
In contrast, M-^jmf/L shows a less tight relation with 
Mdyn (upper right panel), giving M^jmf/L oc A/^^^'^f 
and scatter of 0.03 dex. Comparing the slopes of the two 
relations shows that dark matter/IMF variations con- 
tribute three times as much to the total tilt of the FP as 
do stellar population va riations. 

Figures WM use the iGallazzi et~all (|2005| ) values of 
Mi,jMF/L, but any other choice gives qualitatively sim- 
ilar results: in the FP midplane, lines of constant 
Mi^jMF/L follow lines of constant a, while lines of con- 
stant Mdyn/Mi,jMF follow lines of constant Mdyn- How- 
ever, quantifying the various A-A relations depends on 
the choice of Mi,jMF/L. 

These differences are summarized in Table [TJ The top 
section of the table shows the stellar population and dark 
matter/IMF contributions to the a tilt of the FP, mea- 
sured using the five different estimates of M ^ , tmf /L that 
were discussed in §4.31 The IGallazzi et al.l (j2005| ) values 
are those illustrated in the left panels of Figure [11 values 
using other M^^jmf/L estimates are computed in the 
same way. Notice that, depending on the chosen method 
for estimating M-i,jmf/L, the stellar population contri- 
bution to the a tilt of the FP ranges from 13-70%, with 
our preferred value falling in the middle of the range at 
32%. 

The middle section of Table [T] quantifies the stellar 
population and dark matter/IMF contributions to the 
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Figure 9. The two tilts of ttie Fundamental Plane. Only galaxy bins from the FP midplane are shown. Left: The stellar population 
(upper panel) and dark matter/IMF (bottom panel) contributions to the a tilt of the FP. The stellar population component is a tight 
function of a, with Mi,jMF/L ex o-""^" and scatter of 0.01 dex. The dark matter/IMF component shows more scatter when plotted against 
o", with Miiyn/Mi,j]\[p oc a^'^^ and scatter of 0.05 dex. Comparing the slopes of the two relations indicates that dark matter/IMF effects 
contribute 2/3 of the cr tilt on the FP. Right: The stellar population (upper panel) and dark matter/IMF (lower panel) contributions to the 
total tilt of the FP. Here, the tightness of the relations is the opposite way round: the stellar population term shows 0.03 dex of scatter at 
fixed Mjjy„ with M^jmf/L oc Mjj^^^, while the dark matter/IMF term shows only 0.02 dex of scatter and has M^yn/h[i,ji^ip oc M^„^. 

Dark matter/IMF efl'ects contribute 3/4 of the total tilt of the FP. 

total tilt of the FP, based on the five different estimates of 
Mi,jMF/L. For any given estimate, stellar populations 
contribute less to the total tilt of the FP than they do 
to the a tilt, ranging from 13-58%, with our preferred 
value falling in the middle at 25%. 

For both the a tilt and the total tilt, using the single 
burst values gives the strongest variations in M^^jme/ L. 
This is due to the fact that single burst models likely 
underestimate M^,jmf/L for galaxies with the lowest 
values of M-^jmf/L^ making the trend look stronger 
than it actually is (see section 14. 2p . In contrast, the 
photometry-based deter mina tions of M j, imf/L from 
iBlanton fc RoweisI (|2007[) and lBellet all (|2003l) produce 
very little variation in M^^jmf / L — all galaxy bins have 
nearly identical values — so that stellar population effects 
contribute very little to the a tilt or total tilt using these 
estimates. 

The key point is that, regardless of which method is 
used to estimate M^./mf/^, both stellar population and 
dark matter/IMF effects arc required to explain the tilt 
of the FP. Furthermore, these two components of the FP 
tilt rotate t he plane around d ifferent axes. 

Recently, iTreu et aD (j2010f ) have argued using entirely 



independent methods that Mdyn/J^i-kjMF (which they 
parameterize as an "IMF mismatch parameter" a) varies 
as a function of a^ such that higher-cr galaxies have higher 
a. They find log{Mdyn/M^jMF) oc (1.20 ± 0.25) log cr 
and argue that these variations are either due to differ- 
ences in the IMF as a function of cr, or differences in the 
central profile of the host dark matter halo. The slope 
of the trend they find is somewhat larger than what is 
presented here; their values of M-^jmf are inferred from 
multiband photometry and arc thus best compared to 
the Blanton et al. or Bell et al. values derived here in 
Table [U which give \og{Mdyn/Mi,jMF) oc 0.83 log cr and 
\og{MdynlMi,jMF) oc 0.73 log cr, respectively. 

5.2. The Thickness of the FP 

Having discussed the tilt of the FP in detail by exam- 
ining trends along the FP midplane, we now switch over 
to examining trends through the thickness of the FP. 
Figures [7] and [8] showed that the high- and low-SB slices 
above and below the FP midplane exhibited the same 
quantitative trends as the midplane, but with zeropoint 
shifts in both Mi^jMp/L and Mdyn/Mi,jMF- In this sec- 
tion, we present a detailed exploration of the trends in 
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Table 1 

Stellar Population Contribution to the ^^cr Tilt" , Total Tilt, and the Thickness of the FP 
for Various Methods of Measuring M^,j]\tp 



Gallazzi Kauffmannt Single Burst-f Blanton§ Bell^ 
The a Tilt of the FP 



Alog(A/^,_fMF/^) 
A logo- 

Alog(A/d^„/M^,jMF) 
A logo- 



Fraction from Stellar Pops 



0.30 

0.65 

32% 



0.45 

0.50 

47% 



0.67 



0.12 0.22 



0.29 0.83 0.73 

70% 13% 23% 



The Total (Md^„) Tilt of the FP 



Alog{M^,.iMF/L) 
A log AIay„ 

Alog{May„/AhjMF) 

A log Mdyn 



0.08 



0.24 



Fraction from Stellar Pops 25% 



0.13 


0.18 


0.04 


0.07 


0.19 


0.13 


0.28 


0.25 


41% 


58% 


13% 


22% 



The Thickness of the FP 



Alog(A/^jMF/^) 

Alog/e 

Alog{May„/AhjMF) 
A log I,, 

Fraction from Stellar Pops 



-0.22 

-0.79 

22% 



-0.20 

-0.81 
20% 



-0.54 

-0.48 
53% 



-0.02 -0.02 

-1.00 -0.99 

2% 2% 



These are our 



Note. — 

*Values computed using the lGallazzi et ahl II2005I ) measurements of Af^. /^/F 
preferred values, as discusse d in section 14.31 
fValues computed using the [Kauffmann et al7 l|2003bl ') measurements of Mt,jf,[ p. 
I Values computed using the single burst measurements of Mi, jmf from section [4!T] 
§Values computed using the lBlanton fc RoweisI l|2007l l measurements of M^^jnip. 



^Values computed using the lBell et al.l II2003I ) measurements of M-^j 



both quantities through the thickness of the plane. 

To illustrate this, we adopt the binning strategy illus- 
trated in Figure [H where bins are now defined in a 2D 
cross-section through the FP. Figure [TOl maps Mj^jmf/L 
and Mdyn/Mi,jMF across the FP cross-section. The 
same color scale is used in both panels to facilitate com- 
parison. Mi,jMF/L and Mdyn/Mi,jMF behave similarly 
in this projection; both have maximum values among 
high-cr, low-/e galaxies and minimum values at low a and 
high /e . Lines of constant values run diagonally from up- 
per right to lower left. The main difference between these 
two parameters is the range of variation. The range of 
M*,imf/ L values is much smaller than the range covered 
by Mdyn/M^jMF, particularly in the Ale dimension^ 

We have already quantified changes in M-i,jmf I L and 

^^ This projection of the FP averages over _Re. Fig- 
ures [7HS] showed that, while A/^ imf/L is independent of Re, 
Mdyn/^i*.lMF is not. Thus the Mdy,^/Mi,jMF values mapped 
in Figure [TOl average over a range of M^y^/Mf^jMF values at each 
point in the FP cross-section and are presented for illustrative pur- 
poses only. We have dealt with this variation in all quantitative 



Mdyn/M^^jMF along the plane in Figure IHl Figure [TT] 
makes a similar comparison for the variations through 
the thickness of the plane. To focus on the trends 
through the plane, we normalize out any trends along 
the plane by subtracting mid-plane values of M-^jmf/L 
and Mdyn/Mi,jMF at each bin in (cr, Re). The var- 
ious (cr. Re) bins are represented by gray lines. The 
thick black lines show fits to the ensemble of data values, 
with slopes as indicated. Comparing these slopes shows 
that Mdyn/Mi,jMF varies more than three times as much 
through the thickness of the FP as does Mi,jMF / L. Dark 
matter and/or IMF variations therefore provide the dom- 
inant contribution to the thickness of the FP. 

As with the FP tilt, the relative contributions of the 
stellar population and dark matter/IMF terms depends 
on which estimates of M^^jmf/L are used. These quan- 
titative differences are summarized in the bottom sec- 



presentations of the Mdyn/Mi,j]^.ip-\ogIe A-A relations (Table 
[T] and Figure llll l by explicitly subtracting off the dependence of 

Mclyn/Mi,^IMF on Re. 
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Figure 10. The contribution of (a) stellar population variations and (b) variable dark matter/IMF to the thickness of the FP. Individual 
galaxies are plotted as points in A/g— O" space, as in Figure[2] Color contours indicate Mi,^iMF 1^^ and Mciyn/M*,IMF and are constructed 
in the same way as Figures [THSI The color maps in (a) and (b) are shown on the same scale. Both Mi,jMF /L and M^iyn/Mi^MF vary 
through the thickness of the FP with A/e. The variation in Mdyn/^I*.IMF is much stronger than the variation in Mt,jMp/L and therefore 
dominates the thickness of the FP. 



tion of Table [H Using the iGallazzi et~all (|2005[ ) mea- 
surements of Mi,jMF/L, stellar population effects con- 
tribute only 22% of the tilt of the FP — the remaining 
78% must be due to dark mat ter and/or IMF variations. 
The lKauffmann et al.l ()2003bD measurements give similar 
results. As we saw with the FP tilt, using single burst 
models to measure M-^jmf/L through the thickness of 
the FP results in a larger contribution from stellar popu- 
lation effects (53%) because low values of Mi,jMFlL are 
underestimated, while the photometry-based measure- 
ments lead to very small contributions from stellar popu- 
lation effects (2%) because of the age-metallicity degener- 
acy. Considering the full range of models for Mnji\ip/ L, 
stellar populations contribution 2-53% of the thickness 
of the FP, which means that dark matter/IMF effects 
contribute 47-98%. 

From Table [T] it is clear that, regardless of which 
method is used to determine M-^jmf^ variations in the 
IMF and/or the inner dark matter fraction in galaxies 
are required to explain both the tilt and the thickness of 
the FP. They contribute at least 1/2-1/3 of each effect, 
and it is likely that they dominate both. 

This leads to an important point. To the extent that 
Mi,jMF/L variations through the thickness of the FP 
can be ignored (since t hey constitute only 2 2% of the to- 
tal variation using the IGallazzi et aLl 120051 models), the 
change in /g through the plane must be due to changes 
in Mdyn/Mi,jMF- However, at a fixed point in a and 
Re, Mdyn IS constant by definition. Thus changes in 
Mdyn/Mi,jMF through the plane are primarily due to 
changes in Mi,j]\jp within Re, which means that stel- 
lar mass surface density (S^ja/f) is changing. In other 
words, the low surface brightnesses of galaxies below the 
plane are due to low stellar mass surface densities, not to 
the fact that their stellar populations are dim. Galaxies 
do not move below the plane by fading; they were built 
that way. 
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Figure 11. Top: The change in Mi,jMF /L through the thick- 
ness of the FP. The gray lines represent different bins in {a, 
Re), with the Mt,jMF /L value of the midplane bin subtracted 
off. M^jmf/L decreases by 0.22 dex for every 1 dex increase in 
A/e. Bottom: The change in M^y,-^/Mt,jMF through the thick- 
ness of the FP. Colors and normalization are as in the top panel. 
Mdyn/^-k,lMF decreases by 0.79 dex for every 1 dex increase in 
A/e- Miiy„/Mi,jMp therefore varies more than three times as 
much as AI-^j]\{p/L through the thickness of the FP. 
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Table 2 

Fitting Formulas for Mass-to-Light Ratio Contributions 





As Functions 


of O, Re, 


and Ale* 










log (Mdyn/Lv) 


= 0.85 
(0.01) 


logo-' 


+ 


0.22 
(0.01) 


logR'e - 


1.04 
(0.01) 


A log 4 


+ 


0.62 
(0.01) 


log (Mayn/M^jMF) -- 


= 0.53 
(0.02) 


logo-' 


+ 


0.22 
(0.01) 


logR'e - 


0.80 
(0.02) 


Alog/e 


+ 


0.20 
(0.01) 


log (M^jmf/Lv) 


= 0.31 
(0.01) 


logo-' 


+ 


0.003 
(0.01) 


logR'e - 


0.24 
(0.02) 


Alog/e 


+ 


0.42 
(0.01) 


log (J^.jmf) 


1.49 
(0.01) 


logo-' 


— 


1.20 
(0.01) 


logR'e + 


0.80 
(0.02) 


A log le 


+ 


3.81 
(0.01) 






As Functions of cr, R, 


o , and le 










log (Mdyn/Lv) 


= 2.00 
(0.01) 


logo-' 


- 


0.97 
(0.01) 


iogif4 - 


0.95 
(0.01) 


logl'e 


+ 


0.61 
(0.01) 


log {May„/A'h,iMF) = 


1.43 
(0.03) 


logo' 


— 


0.71 
(0.03) 


logR'e - 


0.74 
(0.02) 


logl'e 


+ 


0.20 
(0.01) 


log i^'hjMF/Lv) 


= 0.57 
(0.03) 


logo' 


- 


0.27 
(0.02) 


logR'e - 


0.22 
(0.02) 


logl'e 


+ 


0.42 
(0.03) 


log (S,,^„j.) 


= 0.60 
(0.02) 


logo' 




0.27 
(0.02) 


logR'e + 


0.74 
(0.02) 


logl'e 


+ 


3.82 
(0.01) 



Note. — o' = ct/150 km s" 
* A log le is defined as A log le 



R'e = Re/2.5 kpc, and /^ ; 
log/e - (1.16 logo- 1.211 



/e/4OOL0 pc-2. 
gRe+ 0.55) (see Figurc[l]l. 



5.3. Fitting Formulas for Mass-to-Light Ratios as 
Functions of a, Re, and A/g 

In this analysis, -we have quantified the variations of 
Mdyn/L, Mdyn/Mi,jMF, and NUjmf/L throughout 3D 
FP space. These results can be summarized by fitting 
each of these parameters as a function of tr, Re, and 
A/g. We assume that each mass parameter is a first- 
order function of logcr, lo gi?e, and A log/ e, then use 
the IDL package mpfit.pro (jMark-wardtll2009l ) implemen- 
tation of Levenberg-Marquardt minimization to find the 
best such solution. Note that because a, Re, and A/g are 
not orthogonal coordinates, the exponents derived from 
explicitly 3D fitting do not in general exactly match those 
presented from ID fits in sections I5.1H5.2I The correct 
quantification of the dependence on each FP parameter 
must be derived from explicitly 3D fits. 

The resulting functio nal forms are given in the top half 
of Tabic H using the iGallazzi et "all (|2005f ) values for 
Mi,jMF/ L. Values in parentheses under each function 
give the \a error in each parameter of the fit, computed 
from the covariancc matrix of the fit. We also include fits 
to the effective stellar mass surface density, defined as 
'^*,IMF = M-tjMp/i'^'^R'i) i'^ analogy "with the effective 
surface brightness. The parameters of the M-i,jmf/L fit 
in Table [2] sho-w that Adi^jMp/L is independent of Re, as 
discussed in section 15.11 For completeness, in the bot- 
tom section of the table we also include fits as explicit 
functions of log le , rather than using the A/g parameter- 
ization. 

Figure [12] illustrates the quality of these fits. In 
each panel, we plot the median value of Mdyn/L, 
Mdyn/Mi,jMF, Mi^jMp/L, or ^i,jMF as measured for 
each bin against the value predicted by the fits in Table 
[21 Panels a-c are plotted on the same scale to facili- 



tate comparison. The standard deviation of the scatter 
about the one-to-one relation is indicated in each panel. 
The fits do an excellent job of reproducing the median 
observed values, with scatter < 0.02 dcx (i.e., < 5% ) 
in all cases. The relations can therefore be used to esti- 
mate values of Mdyn/Mi^jMF, Mi^jMF/L, and Y.i,jMF 
for ensembles of SDSS galaxies from the parameters {a), 
{Re), and {le)- Note that these values are correct in the 
statistical sense only, as the galaxy-to-galaxy variation 
in stellar population properties within a bin cannot be 
estimated from the stacking analysis presented here. 

6. COUPLED DARK MATTER/IMF VARIATIONS AND 
STAR FORMATION HISTORIES THROUGH THE 
THICKNESS OF THE FUNDAMENTAL PLANE 

We have examined in detail the mass-to-light ratios 
of our sample of ^16,000 quiescent galaxies from the 
SDSS. We have identified and quantified the contribution 
of stellar population effects {M^^^mf/L) both along the 
FP and through its thickness. By comparing these mod- 
est variations with the much larger observed variation in 
Mdyn/L, we have argued for substantial additional varia- 
tions in Mdyn/Mi,jMF due to variations in the IMF or in 
the inner dark matter fractions in some galaxies. In this 
discussion, we combine this information with the stellar 
population results of Paper II to evaluate three possible 
physical mechanisms behind the observed Mdyn/Mi,jMF 
variation, with particular focus on variations through the 
thickness of the FP. 

For reference, the correlations between the structural 
properties of galaxies, their mass ratios, and their stellar 
population properties are summarized in Table [31 This 
table combines the results of this work and those of Pa- 
per II. For each of the stellar population properties de- 
rived in Paper II— age, [Fe/H], [Mg/H], and [Mg/Fe]— 
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Figure 12. The quality of fits to (a) M^y^/Lv, (b) M^^jmf / Ly , (c) Mjy„/M^,jMF, and (d) T,i,jMF as functions of cr, Re, and Ale. 
In each panel, the median observed value of the parameter in question is plotted against the predicted value from the fits given in Table [2l 
Solid lines showr the one-to-one relation. The standard deviation about the one-to-one relation is indicated in each panel and is < 0.02 dex 
(i.e., < 5%) for all four relatio ns. The fits given in Table [2] are therefore excellent approximations to the measured values of all four mass 
parameters. These fits use the lGallazzi et ahl I I2005I) values for M^, imf/L. 



Table 3 

Trend Directions for Stellar Population and Structural Properties 



Age 

(Gyr) 



[Fe/H] 
(dex) 



[Mg/H] 
(dex) 



[Mg/Fe] 
(dex) 



{Mq pc"2) 



(Mq/Mq) 



M-^jmf/L Atf 
(Mq/Lq) (Gyr) 



a 


+ 


+ 


+ 


+ 


RJ 


null 


null 


null 


null 


A/e* 


— 


+ 


+ 


— 



+ 



+ 



+ 


+ 


— 


+ 


null 


null 


— 


— 


+ 



Note. — 

* Variations with a along the FP, at fixed Re- 
t Variations with Re along the FP, at fixed a. 
^Variations with Ale through the thickness of the FP, at fixed a and fixed Re 



the table indicates whether that property correlates (+) 
or anti-correlates (— ) with a (along the FP midplane at 
fixed Re), with Re (along the FP midplane at fixed a), 
and with Ale (through the thickness of the FP). The 
table also summarizes the results of this work by indi- 
cating how I]*jj\fF, Mdyn/M^jMF, and IV^jme/L vary 
with a. Re and A/g. Lastly, the table includes an esti- 
mated trend in the duration of star formation {Atsp), 
based on the assumption that star forma tion timescale 
is inversely correlated with [Mg/Fe] (e.g.. iTinslevI 119791 : 



Greggio fc Renzinilll983l:[Worthev et al.1119921: [Matteuccil 



1994HTrager et al.|[2000ai iThomas et al.ll2005[ ). 

The following sections consider three families of mech- 
anisms for producing Mdyn/Mi,jMF variations, focusing 
on variations through the thickness of the FP. For each of 
these, we examine the implications for the stellar popula- 
tion properties and star formation histories of the galax- 
ies, and compare these to the observations. We find that 
the full set of observations is best explained by a model 
that postulates the premature truncation of star forma- 
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Table 4 

The Effects of Variable IMFs 



hh/M^,jMF 
(Mq/Mq) 



A log le 
[Lq pc-2) 



[Mg/Fe] 
(dex) 



[Fe/H] 
(dex) 



[Mg/H] 
(dex) 



Observed 

More very high-mass stars 

More very low-mass stars 



+ 
none 



Note. — 

* Assuming all variation is due to IMF differences, i.e.. 



Mdyn/Mi, = const. 



tion in galaxies with high Mdyn/Mi,jMF (and low /g) — 
those that lie below the FP. 

6.1. Variations in the IMF 

The low-fT galaxies with high A/g (upper left corner in 
FigurefTUb) have log {Mdyn/Mi,jMF) ~ 0, as can be read 
from the color-bar in that figure. This means that the 
stellar mass inferred using a Chabrier IMF can account 
for all of the mass detected dynamically. However, the 
rest of the galaxy bins have log {Mdyn/Mi,jMF) > 0. 
If the entire discrepancy between Mdyn and M-^jnip 
through the thickness of the FP is due to IMF variations, 
then all hut the highest-AIg galaxies must have formed 
with IMFs that produce less light for a given mass of stars 
formed than does the Chabrier IMF. 

There are two ways this could happen: through an ex- 
cess of very low-mass stars that contribute substantially 
to M^, but produce very little light, or through an ex- 
cess of massive stars at early times that create numerous 
compact remnants at late times, again contributing sub- 
stantial mass but little light. 

Variations in the IMF will also change the chemical en- 
richment history of a galaxy and predictions for this can 
be compared to the observed abundances. Table |4] sum- 
marizes the predictions for several different models and 
compares them to the observations. The SDSS galaxies 
show stellar population variations through the thickness 
of the FP such that galaxies with lower A/g have higher 
[Mg/Fe], lower [Mg/H], and lower [Fc/H] than galaxies 
at the same a and Re with higher A/g (Paper II) . These 
trends are shown in the first row of Table HI 

An IMF which produces more massive stars (> 8Mq) 
at early times (and therefore more compact remnants 
at late times), in addition to having higher Mi,/AIi,j]\jF 
will produce a larger number of SNe II relative to SNe 
la. These galaxies should then have higher values of 
[Mg/Fe]. The large number of massive stars will also 
result in higher overall effective yields of both Mg and 
Fe, leading to higher values of [Mg/H] and [Fc/H], as 
summarized in the second row of Table |4l This combina- 
tion of high [Mg/Fe] and higher effective yields docs not 
match the observed stellar population variations through 
the thickness of the FP. 

An IMF with an excess of low-mass stars could also 
produce high values of M^/M^-^mf, such as an IMF with 
Salpeter slope continuing down to low masses rather than 
turning over near IMq, as in the Kroupa and Chabrier 
IMFs. Changing the fraction of stars below ~ ^Mq 
would not affect the relative contributions of SNe II and 
SNe la but it would result in fewer SNe overall, thus sig- 
nificantly reducing the effective yields from both types 
of SNe. These lower yields are consistent with the lower 



values of [Fe/H] and [Mg/H] observed in galaxies with 
high Mdyn/Mi,jMF, but this model does not reproduce 
the increase in [Mg/Fe] seen in such galaxies (see Table 

a. 

In general, it is difficult to simultaneously produce 
higher values of [Mg/Fe] and lower yields in both [Mg/H] 
and [Fe/H] by adjusting the IMF. This is because IMF 
modifications that produce a larger number of massive 
stars will tend to increase both [Mg/Fe] and the effective 
yield. This is not a problem when studying abundance 
trends along the FP, where [Mg/Fe], [Fe/H], and [Mg/H] 
all increase together with increasing a (e.g.. Table [3]), 
but it presents a challenge when interpreting the trends 
through the thickness of the FP. 

A possible way to save the IMF model is to assume 
that a larger proportion of massive stars will result in 
increased levels of SN feedback. Thus, although the 
raw yields will be higher in galaxies with more massive 
stars, a larger fraction of metal-enriched material may 
be removed through intense SN feedback. However, this 
involves an additional assumption about differences in 
feedback between galaxies. In what follows, we assume 
that the observed Mdyn/Mi,jMF variations are domi- 
nated by real differences in dark matter fraction within 

Re- 

6.2. Variations in Inner Dark Matter Fraction: 
Redistributing Stars and Dark Matter through 
Merging 

If dark-to-stellar mass ratio is really varying through 
the thickness of the FP, one of two statements must be 
true. Either galaxies with higher Mdyn/ Mi, are genuinely 
deficient in stellar mass for their halo size, or the over- 
all stellar mass fractions are the same but stars are dis- 
tributed differently such that the dark matter fraction 
within Re is larger We discuss the latter effect in this 
section; the former will be discussed in the next section. 

Substantial dissipation during the galaxy formation 
process will lead to galaxies with centrally concentrated 
baryons, while less dissipation during formation will pro- 
duce more extended stellar distributions. Such a process 
has been proposed by several authors as an explanation 
for the tilt of the FP. Here, we consider whether it might 
also contribute to the thickness of the FP. 

The differential dissipation scenario was proposed more 
than a decade ago to e xplain differences in the kinematic 
struc t ures of galax i es (iBender et al.lll992tlGuzman et al.l 



rT993t ICiotti et al.l Il996t iFaber et al.l 119971 )" More re- 
cently, it has been modeled in the numerical simulations 

^* The dynamical mass estimator M^y^ is predominantly sensi- 
tive to the mass inside Re only. 
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of several groups ( Onorbe et al.ll2005t iKobavashil [20051 : 



iDekel fc Coxl2do6t]Robertson et al.ll2006l ) and invoked as 
an explanation for the tilt of the FP. The amount of dissi- 
pation during galaxy mergers is expected to be largest for 
low-m ass galaxies, due to their higher gas fractions (e.g., 
iKannappan 20041 ). lower surface densities, and relatively 
long timescales for star formation (e.g., iNoeske et al.l 
|2007() . With reasonable assumptions for the depen- 
dence of di ssipation on g alaxy mass, [Robertson et al.l 
(|2006( ) and IDekel fc Co^ (|2006| ) are able to reproduce 
the tilt of the FP as a func tion of M^ through varia- 
tion in Mdyn/Mi,. Recently, iCovingtoiil (|2008| ) has ver- 
ified these results in a more realistic cosmological con- 
text by tracing galaxy properties through semi-analytic 
models of galaxy evolution, while iHopkins et al.l ()2008aD 
have found supporting observational evidence in the 
light profiles of nearby ellipticals. Once the FP has 
been established, further dissipationless merging will 
maintain the FP slope in isolated ellip t ical galaxies 
(iCapelato et al.l Il995t iDantas et al.l 120031: iNipoti et al.l 
120031: iBovlan-Kolchin et al.ll2005HRobertson et al.ll2006t ). 
although this is likely not true for b rightest cluster galax- 
ies (BCGs, e.g.. lLauer erall[2007[ ). 

We have argued above that variations in Af(j,y„/M^ are 
needed not only to reproduce the tilt of the FP, but its 
thickness as well. Suppose that these variations arise 
from different amounts of dissipation in the mergers that 
create different early type galaxies at fixed a and Re . As 
shown in Table |3l the galaxies with lower Mdyn/Mi,jMF 
at fixed (cr, R^) are the galaxies with younger mean ages 
and lower [Mg/Fe], consistent with their having more 
extended star formation histories. In this scenario, dissi- 
pational mergers with their accompanying bursts of star 
formation would have to occur preferentially in galaxies 
at later times, while early type galaxies that had star for- 
mation truncated early would have to have experienced 
a dissipationless merger. 

However, this seems backwards to expectations, as 
galaxies at early times were in general more gas rich. 
Thus earlier mergers might be expected to be more dissi- 
pative, not less. In addition, the trends in Mdyn/M^^jMF 
and in stellar population properties continue unchanged 
down to small galaxies with a « 100 km s^^ (see Figures 
ITHTUI and Figures 7-10 of Paper II), where SO galaxies 
are prevalent. Their low values of Mdyn/M^-jMF require 
them to have experienced a substantially dissipational 
merger under this scenario, after which they would have 
to regrow a disk and further quench star formation in 
the residual disk. It is not clear that this merger-driven 
mechanism can produce these smooth trends down to the 
modest galaxy masses included in this sample. 

6.3. Variations in Inner Dark Matter Fraction: Low 
Efficiency Star Formation 

The second scenario for variations in dark matter 
fraction consists of global differences in the stellar 
and dark matter mass fractions within a halo. In 
this scenario, a significant fraction of the baryonic 
material in some haloes does not get converted into 
stars. We will refer to this as a low "conversion ef- 
ficiency" for the baryons. This conversion efficiency 
is known to vary with mass, such that Mi,/Mhaio has 
a maxi mum around the characteristic ga laxy luminos- 
ity (L*, iBenson et all 120001: iMarinoni fc Hudson, ,2002; 



Ivan den Bosch et al.ll2003HZaritskv et al.ll2006f ) and from 
there decreases toward lower masses (e.g., dwarf galax- 
ies) and toward higher masses (e.g., BCGs and galaxy 
clusters) . 

The galaxies in our sample are on the higher-mass 
side of this curve and should therefore have Mhaio/M^, 
increasing with mass, i.e., their conversion efficiency of 
baryons into stars should decrease with increasing mass. 
The physical mechanism driving this change is thought 
to involve the truncation of star formation at some point 
in the evolution of the central galaxy, such that the 
galaxy halo continues to accrete mass but the incom- 
ing baryons are prevented from forming stars. It is very 
interesting that Mdyn/Mi^jMF is observed to increase 
with Mdyn along the FP (e.g.. Table |3]). This increase 
in Mdyn/Mi^jMF could well be the same phenomenon as 
the overall increase in Mhaio/Mi,, both of which appear 
to scale directly with mass. This suggests that the con- 
version efficiency within i?e may mimic the conversion 
efficiency of the parent halo as a whole. 

As we did in the previous section with dissipational 
merging, we can invoke a process thought to cause 
Mdyn/Mi,jMF variation along the FP to try to explain 
the variations through the thickness of the FP as well. 
The time at which star formation is truncated may vary 
among galaxies at the same a and Re, such that some 
galaxies have their star formation shut down prema- 
turely, while other galaxies continue forming stars, accu- 
mulating more stellar mass, and further enriching their 
ISMs for extended periods of time. 

This scenario of variable truncation times finally pro- 
duces the consistent picture we have been looking for to 
explain all the trends observed through the thickness of 
the FP. Galaxies at fixed cr and Re (i.e., fixed Mdyn) 
which are prematurely truncated wind up with shorter 
timescales for star formation (i.e., higher [Mg/Fe]), older 
mean ages, lower total stars formed (lower M^, and lower 
^*.imf), and incomplete processing of metals (lower 
[Fe/H] and lower [Mg/H]). In contrast, galaxies that 
are not truncated continue forming stars longer, result- 
ing in longer star formation timescales (lower [Mg/Fc]), 
younger mean ages, more total stars (higher Af^ and 
higher E^/mf), and more complete metal enrichment 
(higher [Fe/H] and [Mg/H]). 

One possible truncation mechanism for satellites is 
the rapid quenching of star formation when a satel- 
lite f alls into a massive halo and is stripped of gas 
(e^ lGunn fc GottI 11971 ILea fc De Yml^Wi^ IGisleii 
Il976f ) . This explanation predicts that galaxies with high 
Mdyn/Mi,jMF (and therefore high [Mg/Fe]) at a given a 
would tend to be satellites in more massive haloes today. 
However, there is some evidence that this typ e of quench- 
ing p roceeds slowly in massive galaxies fe.g.. lWolf et al.l 
120091 ) and may not be adequately abrupt to produce 
the high [Mg/Fe] seen in massive galaxies with high 
Mdyn/Mi,jMF- Moreover, many of our sample galaxies 
are central to their haloes — 90.1% of our sample galaxies 
are id entified as the mos t massive galaxies in their haloes 
in the I Yang et al.l ()2007f ) group catalogs — and thus satel- 
lite quenching processes are not relevant for them. 

Another possible truncation mechanism is very pow- 
erful feedback, in which a substantial quantity of 
the interstellar medium (ISM) is heated and removed 
from the galaxy potential well before it can form 
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stars. Both starburst-drivcn outflows (with veloci- 
ties up to ^ 600 km s~^ iHeckman et alj 120001) and 
AGN-driven outflows (which can reach velocit ies in 
excess of ~ 1000 km s~^ I Trump et al.l I2006D may 
be able to unbind gas from galaxies. Theoreti- 
cal frameworks for AGN feedback in p articular have 
been provided by many autho rs fe.g.. | Granato et al.l 
2003; 'Scanna pieco fc Obi [20p: iDi Matteo et al.l 12005^ 
Springcl et al. 120051: iHopkins et al.ll2008bD . Real hfe ex- 
amples of AGN feedback could include the two massiv e 
z w 3.5 galaxiaes discussed in iNesvadba et al.l (J2007D . 
They present optical and radio observations the strong 
radio jets observed in these systems and calculate that 
the mass-loading of the jet material may represent sev- 
eral times 10^'^ Mq of ISM. Such objects may be evidence 
of powerful AGN feedback in action. SN-drivcn out- 
flows may also be able to quench star formation in early 
type galaxies (e.g.. lPipino et al.ll2008l ). although there is 
some evidence that the outflow velocities in identifiable 
post-starburst galaxies require the presence of an AGN 
iTrcmonti ct al. 200j). 

A variation on feedback-induced truncation could 
be supernova feedback in smaller galaxies that then 
assemble into larger galaxies. Galaxies with high 
Mdyn/M-i,jMF could be galaxies that did not form the 
majority of their stars in situ, but instead assembled hi- 
erarchically at early times from many smaller galaxies. 
These small galaxies might have shallow enough potential 
wells that supernova f eedback can play an huportant role 
in ren ioving baryons (iWhite fc Reeslll978| jDekel fc Sillj 
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[19861 IWhite fc Fre'^ 119911: IBenson et al.il2003D . This 
process could shut down star formation at earlier times 
and lead to the observed older ages and higher [Mg/Fe] 
of their massive descendents, as compared to galaxies at 
the same a with lower Mdyn/J^*,iMF, which may have 
formed more of their stars in situ in massive haloes. A 
galaxy which formed in many smaller pieces would also 
have lower metallicity, both because its ISM experienced 
fewer generations of stellar processing resulting in a lower 
effective yield, and because SN feedback would be more 
effective at removing metal-enriched gas from the galaxy. 
In this scenario, the galaxies that form in smaller pieces 
would also have to form at earlier times (to match the 
observed age trends). 

Yet another mechanism is que n ching by massive haloes 
(e.g., iBirnboim fc Dekeli l200l iCattaneo etHI l2008h . 
when halos pass over a critical mass threshold {Merit) 
and shock-heated gas accreting onto them can no longer 
cool efficiently. Variations through the FP could be due 
to stochastic variations in the halo mass assembly history, 
if galaxies with the same a and Re today exist in halos 
that passed over Merit at different times, or if mass was 
accreted differently onto these halos, affecting the cooling 
rate. This would mean that halo mass is not perfectly 
coupled to the structure of the visible galaxy. Also, halo- 
quenching scenarios need to be coupled with some kind 
of maintenance mode to suppress star formation, such as 
"radio- mode" AGN activity (e.g.. lCroton et al.ll2006h . 

All of these scenarios of truncated star formation 
are qualitatively consistent with the higher values of 
[Mg/Fe], lower effective yields, and hypothesized shorter- 
duration star formation histories of high Mdyn/Mi,j]\jF 
galaxies. Indeed, low-efHciency star formation due to 



truncation is an attractive solution to simultaneously 
achieving high [Mg/Fe], low effective yields, and lower 
total production of stars. Semi-analytic models of galaxy 
formation may provide clues as to which, if any, of these 
mechanisms can produce quantitative agreement with 
the observed galaxy properties. 

7. CONCLUSIONS 

In Paper II, we demonstrated that the stellar popu- 
lation properties, and therefore the star formation his- 
tories, of quiescent early type galaxies form a two- 
parameter family. These map onto a cross-section 
through the Fundamental Plane. In this paper, we have 
explored the associated mass-to-light variations in 3D FP 
space. This has allowed us to map out the various com- 
ponents of mass-to-light variation both along the FP and 
through its thickness, with the following results: 

1. We map Mdyn/L, Mdyn/M^jMF, and M^jmf/L 
throughout 3D Fundamental Plane space. These 
quantities can be well-approximated by first-order 
functions (i.e., hyperplane fits) of logcr, logi?e, and 
Alog/g. These fits are provided in Table [2l 

2. We confirm that the variation in M^^jmf/L due 
to stellar population variations cannot account for 
the tilt of the FP, as demonstrated by many pre- 
vious authors. Variations in the dark matter frac- 
tion within Re (i.e., Mdyn/M-t) or in the IMF (i.e., 
Mi,/Mi,jMF) are required to explai n the tilt of the 
FP. U sing the M^,jmf estimates of iGallazzi et al.l 
()2005l ). stellar population effects contribute only 
1/3 of the tilt of the FP, while dark matter/IMF 
effects contribute the remaining 2/3. 

3. The stellar population and the dark mattcr/IVIF 
contributions to the FP tilt rotate the plane around 
different axes in the 3D space, with M-i,j]\[f/L oc 
f{(7), while Mdyn/M^jMF oc g{Mdyn)- 

4. Although Paper II showed that stellar populations 
vary through the thickne s s of the FP (in agree ment 
with iForbc s ct al.J_ |l99l IWuvts et all 12004 and 
iGargiuloet al. 2003), the associated variations in 
Mi,j]\{F/L alone cannot account for the thickness 
of the FP. Variations in the dark matter contribu- 
tion within Re (i.e., Mdyn/Mi,) or in the IMF (i.e., 
Mi,/Mi,jMF) are required to explain the thickness 
of the FP (a s well as the ti l t). U sing the Mi,jMF 
estimates of IGallazzi et al.l (j2005f ). stellar popula- 
tion effects contribute only 22% of the thickness of 
the FP, while dark matter/IMF effects contribute 
the remaining 78%. 

5. Because stellar population effects contribute only 
~ 22% of the thickness of the FP, the observed 
Mdyn/L variations through the plane are domi- 
nated by variations in T,i,jMF- This means that 
galaxies do not move below the plane by fading; 
they either form initially with lower E^,jj\./f or 
move there due to structural changes in the stel- 
lar mass distribution. 

6. Combining the results of this work with the re- 
sults of Paper II, wc find that variations in the 
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dark matter/IMF contribution through the thick- 
ness of the FP are associated with differences in 
galaxy star formation histories. The correlation is 
such that galaxies with higher Mdyn/Mi,jMF also 
have higher [Mg/Fc], older ages, and lower mctallic- 
ities than their lower- Mdyn/ M^^j m f counterparts 
at the same a, which may indicate that the higher- 
Mdyn/Mi^jMF galaxies have experienced shorter 
duration star formation. Possible physical mech- 
anisms for producing h[gh-Mdyn/Mi,jMF and the 
associated stellar population variations include: 

a) A top-heavy IMF in some galaxies, which pro- 
duces more SNe II relative to SNe la and 
therefore higher [Mg/Fe]. However, this sce- 
nario does not simultaneously produce the 
observed lower [Fe/H] and [Mg/H] in these 
galaxies. Invoking IMF variations to ex- 
plain the observed Mdyn/J^*jMF variations 
would require an additional physical mech- 
anism, such as enhanced SN feedback that 
removes a substantial fraction of the metal- 
enriched material. 

b) Higher dark matter fraction within i?e due to 
low-dissipation merging in galaxies that have 
already quenched their star formation. This 
scenario is plausible, but the predicted trend 
with stellar population age seems backwards 
to expectations since early-forming galaxies 
are naturally more gas-rich, not less. 

c) Low conversion efficiency of baryons into stars 
through the premature truncation of star for- 
mation. There are various possible mech- 
anisms for this trunction, including ram- 
pressure stripping in satellite galaxies, power- 
ful AGN feedback, massive halo quenching, or 
low-mass progenitor galaxies that are strongly 
affected by supernova feedback (as compared 
to galaxies that form most of their stars in 
situ in deep potential wells). This trunca- 
tion scenario is appealing because it simulate- 
nously matches the lower stellar masses, older 
ages, higher [Mg/Fe], lower [Fe/H], and lower 
[Mg/H] observed in the h\g}\- Mdyn I M^^jmf 
galaxies at fixed a. 

7. There is an intriguing parallel between the known 
increase of M halo /Mi, in large haloes and the in- 
crease in Mdyn/Mi^jMF along the FP observed in 
this galaxy sample. This suggests that the con- 
version efficiency of baryons into stars inside Re 
mimics the overall conversion efficiency for the par- 
ent halo. The truncation of star formation in mas- 
sive galaxies is widely accepted as the explanation 
for the increase in Mhaio/Mi, with galaxy mass, al- 
though the truncation mechanism is debated. It 
may be that this same truncation mechanism also 
produces differences in conversion efficiency inside 
Re, and that this effect drives the observed varia- 
tion in Mdyn/Mi,jMF through the thickness of the 
FP as well. 

These results have focused on quiescent galaxies only, 
excluding early type galaxies with ongoing star for- 
mation, Scyfcrt, or LINER activity. Seyfert and 



LINER hosts are known to have systematically younger 
stellar p opulation ages than th e ir quiescent coun- 
terparts (iKauff mann et all l2003al: iGraves et all 120071: 
iSchawinski et a l. 2007), but further analysis is required 
to determine whether they follow the other stellar popu- 
lation and AI/L trends explored here for quiescent galax- 
ies. 

Taken as an ensemble, the results presented in this se- 
ries of papers illustrate a connection between the star 
formation histories and the present-day mass structures 
of early type galaxies. As such, they represent a strin- 
gent test for models of galaxy formation and a challenge 
for the next generation of semi-analytic models. They 
also highlight the emerging importance of abundances in 
testing model star formation histories. Paper IV makes 
extensive use of this abundance information to expand on 
and quantify the premature truncation model proposed 
here. 
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Figure 13. Three examples of weak [Olll]A5007 omission in the stacked spectra. In each panel, the thick black line shows the data, the 
thin gray line shows the BC03 model with the closest match in age and [Fe/H], and the lower dark gray line shows the difference between 
the data and the model. The vertical line indicates the location of the [Olll]A5007 emission line, which is weak but clearly detected in all 
cases. The data have, from top to bottom, logcr = 2.06, logcr = 2.14, and logcr = 2.23, with the models smoothed to match the combined 
SDSS spectral resolution and intrinsic velocity dispersion of the data. The BC03 models have a fixed abundance pattern, while many of 
the sample galaxies have substantially super-solar [Mg/Fe] (as labelled). This causes a mismatch between the data and models around the 
5190 A Mg absorption feature but does not seriously bias the fit to most of the spectrum, which is dominated by Fe absorption features. 
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APPENDIX 
CORRECTING FOR LOW-LEVEL H/3 EMISSION INFILL 

The galaxies presented here are selected to contain no ionized gas emission by requiring that the strong emission 
lines at Ha and [Oii]A3727 be undetected at the 2cr level. This results in a sample of quiescent galaxies that lack both 
significant ongoing star formation and AGN/LINER activity. 

The individual SDSS spectra have relatively low S/N. Our sample therefore contains galaxies with weak emission 
below the detection threshold in individual galaxies. Weak H^ emission will partially fill in the H/3 absorption feature 
used to determine galaxy ages, causing the stellar population modelling to return ages that are too old. Even small 
errors in H/3 can result in significant overestimates of the galaxy ages; at 10 Gyr, 0.1 A in H/3 corresponds to an 
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Figure 14. (a) The EW of [Olll]A5007 as a function of the offset in [Fe/H] between the measured [Fc/H] of the data and the closest- 
matching BC03 model. Mismatches in [Fc/H] between the data and the model cause biases in the estimates of [Olll]. We use a third-degree 
polynomial fit (solid line) to fit the trend in [Olll] with A [Fc/H] and correct the [Olll] measurements to the values that they would have if 
A [Fe/H] = for the entire sample (b). The mean and standard deviation of the corrected [Olll] EW measurements are given. 



age difference of roughly 1 Gyr ()SchiavoDll2007D . However, if weak emission can be detected in the high-S/N stacked 
spectra, a correction for this effect can be applied to the measured values of H/3, resulting in more accurate age 
measurements. 
We do this by measuring the equivalent width (EW) of emission in the [Oiii]A5007 line and using the relation 

AH/3 = 0.7EW([Oiii]A5007) (Al) 

suggested bv lTrager at al.l (|2000bD . Using the sample of I Graves et al.l (|2007D . we have confirmed that this relation is 
a good approx imation for the l ine ratios of LINERs, which likely make up the majority of weak emission-line objects 
in our sample (|Yan et al.ll2006l see also Figure 7b of that work). 

To measure EW[Oiii], we restack all of the composite spectra without smoothing the spectra to cr = 300 km s~^. A 
weak, narrow emission feature should be easier to detect if it is not smoothed up to the maximum velocity dispersion 
of the sample. In practice, the width of [Oiii] appears to track the stellar a closely, such that the emission in high-cr 
galaxies is substantially broader than the emission in low-cr galaxies, although the integrated EW is similar. 

For each stacked spectrum, we identify the BC03 model with the closest match to the age and [Fe/H] values measured 
in our stellar population analysis. The BC03 models sample parameter space fairly densely in age, but only sparsely 
in [Fe/H], so that the age mismatch between data and model is less than 0.02 dex, while the mismatch in [Fe/H] can 
be as much as 0.2 dex. Three examples with different a values are shown in Figure [Ql We divide out the continua 
of the data and the model, smooth the model to match the combined SDSS resolution and intrinsic a of the stacked 
spectrum, and compare the data (red line) and model (black line) in the 4800-5200 A range. 

The models are an excellent match to the data except in two locations: [Oiii] emission, though weak, is clearly visible 
at 5007 A in the data but is not included in the models, and the data show stronger absorption in the IVIg feature at 
5190 A than is predicted in the models. The latter effect is due to the fact that the BC03 models have a fixed solar 
abundance pattern, while early type galaxies typically have super-solar []VIg/Fe]. The larger the measured value of 
[]VIg/Fe], the bigger the mismatch between data and model at 5190 A. The mismatch in the abundance pattern does 
not seriously bias the fit to the rest of the spectrum, which is dominated by Fe absorption features 

To measure [Oiii], we compute the difference spectrum between the da ta and the model (blue line) and integrate 
the residual [Oiii] flux, using the bandpass and continuum definitions of lYan et al.l ()2006l ). The detection of [Oiii] 
emission implies that il/3 emission is also present, filling in the stellar H/3 absorption line. This H/3 infill is not visible 
because the chosen BC03 model has been matched to the mean age of the uncorrected stacked spectrum. 

As mentioned above, the BC03 models are only computed for a limited set of [Fe/H] values and are therefore not 
always well-matched to the data. JVIetallicity offsets between the data and the model can bias the measurement of 
EW([Oiii]). Figure [Tib shows the measured value of EW([Oiii]) as a function of the offset between [Fe/H] as measured 
in the data and the closest-matched BC03 model. Compared to galaxies that are well- matched to the models (A [Fc/H 
« 0), galaxies for which the best- fitting model is at higher metallicity (A [Fe/H] < 0) or lower metallicity (A [Fe/H 
> 0) produce EW([Oiii]) measurements that are significantly overestimated or underestimated, respectively. 

To correct our EW([Oiii]) measurements to the value they would have if a perfectly-matched model were available, we 
fit a third-degree polynomial to the data in Figure fT4k and subtract out the trend with A [Fe/H] to produce corrected 
EW([Oiii]) values (Figure fT4b). The corrected EW([Oiii]) measurements have a roughly Gaussian distribution with 
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Figure 15. (a) A comparison of the g — r colors derived from stellar population modelling with the observed g — r colors as measured 
from photometry matched to the 3" SDSS spectral fiber aperture. Model g — r colors are derived by matching the age and [Fe/H] measured 
for the stacked galaxy spectra to BC03 models, as illustrated in Figure |^. The line shows a one-to-one relation, with shading indicating 
the typical la scatter expected due to absorption line measurement errors. Data with substantially larger errors are shown as open circles. 
For high-(T galaxies, the model colors are an excellent match to the observed colors, with observed colors showing < 0.02 mag of internal 
reddening that is not included in the models. The discrepancy is larger for low-c galaxies, suggesting that single burst models systematically 
underestimate the mean stellar population ages for younger galaxies (see text for details), (b) The same as (a), but with models derived by 
matching the age and [Z/H] measured in the stacked spectra. Matching [Z/H] instead of [Fe/H] results in model colors that are substantially 
too red. 

mean /i = 0.22 and standard deviation a — 0.07. There is no correlation between the corrected EW([Olll]) and galaxy 
a, Re, le, A/e, age, [Fe/H], [IVIg/Fe], or mismatch between the measured age and the closest B C03 model age. 

We use the corrected EW([Oiii]) measurements to estimate a correction to H/3 using cquation lAll These corrections 
are subtracted from the measured H/3 line strengths. We then rerun the stellar population analysis on the corrected 
data. Correcting for emission infill produces ages that are ^ 0.12 dex lower than the uncorrected values, while [Fc/H] 
values increase by ~ 0.06 dex. 

These corrected ages are of course different from the original ages used to fit the BC03 models and measure [Oiii] 
emission. To check that the EW([Oiii]) values are reasonable, we iterate the process, using the corrected ages to fit 
BC03 models to the stacked spectra and rcmeasure EW([Oiii]). The resulting values of EW([Oiii]) are typically 0.02 A 
smaller (with standard deviation 0.02 A) than the original estimate, which translates into H/3 differences of < 0.015. 
These differences are negligible. 

COMPARING OBSERVED AND MODELED GALAXY COLORS 

By comparing the color estimates from Figure |4}d to the observed galaxy colors, we can check that this method 
of matching BC03 models to the data gives reasonable estimates of the stellar population properties of the sample 
galaxies. The spectroscopic data sample only the region of the galaxy covered by the SDSS spectral fibers, but many 
galaxies show significant color gradients. We therefore compare the predicted model g~r colors from the spectroscopic 
data to the observed galaxy colors as measured from SDSS photometry matched to the 3" SDSS spectral fiber aperture. 
The fiber magnitudes arc downloaded from the SDSS Catalog Archive ServeiH and K-corrccted to z = 0. We use 
the median value of the g — r fiber color from all the constituent galaxies in the corresponding bin to determine the 
"observed" g ~ r fiber color of each stacked spectrum. 

Figure 115b shows the comparison between the observed g — r fiber colors and those derived from the spectroscopic 
single burst models using the curves in Figure IHd. The various colors indicate the values of a in each bin, as labelled. 
The solid line shows a one-to-one relation between observed and modeled colors, with the shaded region showing the 
typical ItT error on the model colors as propagated through the stellar population modeling. Bins whose stacked 
spectra have considerably lower S/N (i.e., errors in the model colors more than twice as large) than the typical bins 
are shown as open circles. For the high-cr galaxies, the model colors are an excellent match to the observed colors. The 
observed colors are offset not more than 0.02 dex redder than the model colors derived for these galaxies, which may 
indicate low-level internal dust reddening within the sample galaxies. No reddening has been applied to the models. 

At lower values of cr, the discrepancy between observed and modeled colors increases; the galaxies are observed 
to be 0.03-0.04 dex redder than the model colors. This systematic trend, wherein the models predict colors for the 
bluest galaxies that are systematically too blue, highlights the probable bias inherent in using single burst models. 
The Balmer absorption lines are highly sensitive to a small "frosting" of young stars — more so than the total galaxies 
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colors-and will therefore give age estimates that arc systematically younger than models with composite star formation 
histories. These will then lead to predicted colors that arc too blue and, as we show in section|4l to A/^,/ji/i?/L values 
that are too low. 

Given these considerations, the colors derived from the modelling illustrated in Figure |4] appear to do an excellent 
job of matching the observed galaxy colors. Furthermore, they reveal the effect of an expected systematic bias in the 
modeling process. 

A further important issue is also clarified by matching the model and observed colors. The BC03 models are 
constrained to use the solar abundance ratios at all metallicities. However, ma ssive early type gala xies are known 
to h ave super-sola,r enhancements of a nu mber of elements, including Mg (e.g., iWorthev et al.|[l992l ) and C and N 
(e.g.. lGraves et al1l2007t [Smith et al.ll2009l ). Since a mismatch between the galaxy abundance patterns and the model 
abundance patterns is therefore inevitable, it is not obvious a priori whether one should attempt to match [Fe/H] 
between the data and the models, or whether it would be more appropriate to match total mctallicity ([Z/H]). 

Figure [T5k was constructed by matching [Fe/H] between the data and the models. Figure [TSb shows results when 
the match is instead based on [Z/H]. For the observed galaxies, [Z/H] is calculated based on the assumption that the 
a-elements O, Ne, Na, Si, S, and Ti all track Mg, that the Fe-peak elements Cr, Mn, Co, Ni, Cu, and Zn all track 
Fe, and that Mg, Fe, C, N, and Ca have the a bundances measured using EZ_Ages. All other elements are assumed to 
have solar abundances (see [Graves et al.]|2007l section 5.3). 

Matching models to the data on the basis of [Z/H] clearly does not do as well as matching based on [Fe/H]. The 
modeled colors are substantially redder than the observed colors for the high-a galaxies. 

The good match achieved between the observed a nd m odeled colors using [Fe/H] motivates us to use the same 
models to measure Mi^jMp/L. We do this in section |4J1 keeping in mind the caveat that single burst models tend 
to produce colors that are too blue for the bluest galaxies and will similarly underestimate Mi,jMF / L for the lowest 
Mi,jMF/L galaxies. 



